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PREPACK 


Of all the celestial bodies whose motions have formed the 
subject of the mvestigations of astronomers, the Moon has 
always been regarded as that which presents the greatest 
diflSculties, on account of the number of inequalities to which 
it IS subject, but the frequent and important apphcations of 
the results render the Lunar Problem one of the highest in- 
terest, and we find that it has occupied the attention of tlie 
most celebrated astronomers from the earliest times. 

Newton’s discovery of Universal Gravitation, suggested, it 
IS supposed, by a rough consideration of the motions of the 
moon, led him naturally to examine its application to a moie 
sevoie explanation of her disturbances , and his Eleventh Section 
IS the first attempt at a theoretical investigation of the Lunar 
inequalities The results he obtained were found to agiee 
very nearly with those determined by observation, and afforded 
a remarkable confirmation of the truth of his great principle , 
but the geometrical methods which he had adopted seem 
inadequate to so complicated a theory, and lecoursc has been 
had to analysis for a complete determination of the disturbances, 
and for a knowledge of the tiuc oibit 
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The following pages will, it is hoped, form a propei lutio- 
duction to moie recondite works on the subject the difficulties 
which a pel son entering upon this study is most likely to 
stumble at, have been dwelt upon at considerable length, and 
though different methods of investigation have been employed 
by different astronomeis, the difficulties met with are ncaily 
the same, and the principle of successive approximation is 
common to all In the present woik, the approximation is 
carried to tlie second order of small quantities, and this, 
though far from giving accuiate values, is amply sufficient for 
the elucidation of the method 

The differences m the analytical solutions arise fiom the 
various ways m which the position of the moon may be 
indicated by alteimg the system of cooidinates to which it 
IS referred, or again, in the same system, by choosing diffeient 
quantities for independent variables 

D’Alembert and Clairaut chose foi coordinates the pio- 
jcction of the radius vector on the plane of the ecliptic and 
the longitude of this projection To fonu the diffeiential 
equations, the tme longitude was taken foi independent 
variable 

To determme the latitude, they, by analogy to Newton’s 
method, employed the differential vaiiations of the motion of 
the node and of the inclination of the orbit. 

Laplace, Damoiseau, Plana, and also Heischel and Airy 
m then more eleinentaiy works, have found it more con- 
venient to express the variations of the latitude directly, 
by an equation of the same foim as that of the ladius 
vectoi 
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Lubboik cind Pontecoiilantj taking the same cooidinates of 
the moon s positionj make the time the independent variable , 
and wlien it is desiicd to caiTy the approximation to a high 
oidci, this method ofibis the advantage of not requiring any 
revcision of sciies 

I oisboii proposed the method used m the planetary theoiy, 
/ that IS, to doteimiiie the vaiiation m the elements of the 
moon s 01 bit, and thence to conclude the corresponding varia- 
tions ot the radius vector, the longitude, and the latitude. 

llio selection of the method followed in the present work, 
which IS the same as that of Airy, Hcrschel, &c , was made 
on account of its simplicity moreover, it is the method which 
has ohtamod in this umveisity, and it is hoped that it may 
pio\c ot seivicc to the student in lus reading foi the exami- 
nation toi Ilonouis In fuitheiance of this object, one of the 
chapters (the sixth) contains the physical intcipietation of the 
vanoiis impoitaut teims in the ladius vectoi, latitude, and 
longitude ^ 

Tlic seventh cliaptcr, oi Appendix, contains some of the 
most interesting icsults m the terms of the higher ordeis, 
among winch will be found the values of c and g completely 
obtained to the thiid oidei 

^J^'lie last chapter is a biief histoiical sketch of the Lunar 
Ihobicm up to the time of Newton, containing an account 
of the discoveries of the several inequalities and of the methods 
by which they were represented, those only being mentioned 
which, ,is tlic thooiy has since veiified, were leal onward 


^ Sf’o the Hopoit oi the “ lioiul of MathcmutiCiil Studies” for 1850 
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HtopB. Tlio pcnisal of this chapter will shew to what extent 
wo aio mdebtod to our gicat philosophci ; at the same time 
we cannot fail being impiessed with revcience for the genius 
uiul pciseverance of the men who preceded him, and whose 
idaboratc and multiplied hypotheses wei’e in some mcasuie 
noeessaiy to the discovery of his simple and single law. 

I take this opportunity of acknowledging my obligations 
to ftovcral friends, whose valuable suggestions have added to 
t!ic utility of the woik 


HUGH GODFRAY 


lO^A, 1853 
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LUNAR THEORY. 


CHAPl'EK I. 


INTRODUCTION 

Before proceeding to the consideration of the moon’s 
motion, it mil be desiiable to say a few woids on the law of 
attractions, and on the peculiar circumstances which enable us 
to simplify the present mvestigation 

1 The law of universal gravitation, as laid down by 
Newton, is that “ Emery pat tide in the universe attracts every 
other particle, with a force varying directly as the mass of the 
attracting particle and inversely as the square of the distance 
between them ” 

The truth of this law cannot be established by abstract 
reasoning ; but as it is found that the motions of the heavenly 
bodies, calculated on the assumption of its truth, agree more 
and more closely with the observed motions as our calculations 
are more strictly performed, we have every reason to consider 
the law as an established truth, and to attribute any slight 
discrepancy between the results of calculation and observation 
to instrumental errors, to an incomplete analysis, or to our 
Ignorance of the existence of some of the disturbing causes. 

Of the last cause of deviation there is a remaikable instance 
in the recent discovery of the planet Neptune, for our acquaint- 
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ance with which as one of tlic bodies of oiii we are 

indebted to the pcituibutioiis it loodueed iii the taUuiatiul 
places of the pluiet ITuums Tlic‘se piaiuihations too 

great to be attiibuted wholly to eiiois of instrmoeiits ut ot 
calculation, and tliciefoie, eith<‘i tlie low of lutlvn^al 
tatioii was licie at fault, oi som<‘ unknown bodv wa^ tliNtinhnie; 
the path of the planet 'riiis last snpiiosition, in the pow<ulul 

hands of Messis Adams and la‘ Vniui, hsl to the deteetion 

of Neptune by solving the ddlieult niMtst* problem, \i/, 
Given tlie pcrtiubations enuseil by bod\ , detennnn , on the 
assumption oi the tiuth of New^toiTs law, tin* iubit and position 
of the distuibing bcKly. 

Evidence so strong as tins foues ns to admit the eoueit- 
ness of tlie assumption, and 'we must now timsuhr how this 

law, coinhmcd with the laws of motbm, will niabl< ns to m- 

vostigatc tlie ciieiimHtaiK*(‘H of tin* nutoids motuiu, and to u^sign 
hci position at any turn* wlum <,jbs*»nation h.is tmmdusl tlie 
leqiuHitc data. 

2 . The proldein in its pu sent fbnn wtmld be t»ne nf 
extreme, If not msunnoiintable diHitult\ it we bat! to take 
into acoomit simuUaneou^y tlie uetums ft! the eaitb, sun, 
planets, &e. on the moon; hut futunately tin* eaitirs altim- 
tion, on account of its proMinitN, is mmb gitatn than the 
distil) l)iuq] foiee of the sun or of any planet, tin sc disturb- 
ing forces being so small c*ompared with tlu^ ahsolute foMc of 
the earth, that the Hquiirew and products of tlie eflMtH tliey 
produce may be neglected, eveept in extienu’ la^es* and tliiie 
18 a principle, called the ‘^pinietjib* of the sup» rpootinn of hiu ill 
motions,” whieli shows that in such a case the diUmbmg foucH 
may he consideKsl soparatidv, and the algebraic sum of the 


^ It had been wen b) l>i. Lnnumt n! Munnb, out Imm *|h 

known to bo urUuKt, **holui SjMun, by J U Jhad/* 

t Smoe Ibe sun ntti.ulh bnib th** * itth nud laonui 0 ? fb o itiO Oh 
effects on tho moc,m*H inotum rcb»lo#lj Oi iJn i ai^fb or fln^ ticUmbinv? bau 
will not be tlu Mime an the ibsolute ban’ »m < Uhtr I ho wili hr biltv 

mvestigatod m Arts, (0) and fid), 
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distmbances so obtained will be the same as the disturbance 
due to tlie simultaneous action of all the forces 


Principle of Superposition of Small Motions 

3 Let a particle be moving undei the action of any number 
offerees some of which are very small, and let A (fig 1) be the 
position of the particle at any instant Let two of these small 
forces be omitted, and suppose the path of the particle 

undci the action of the lemaming forces to be AP in any 
given time 

Let AP^ be the path which would have been described in 
the same time if also had acted, AP^^ differing very shghtly 
from AP^ and PP^ bemg the disturbance 

Similarly, if mstead of w, had acted, suppose AP^ to 
have represented the disturbed path {AP, AP„ AP^ are not 
necessaiily in the same plane, nor even plane cuives), PP 
bemg the disturbance * 

Lastly, let AQ he the actual path of the body when both 
m, and act Join P,Q 

Now, since the path AP^ very nearly comcides with AuP, 
the disturbance P,Q, due to the action of m,, on the path AP„ 
can diffei in magnitude and direction from the disturbance PP^, 
due to the action of the same force on the path AP, only by 
a quantity of the first order compared with or of the 
second order compared with AP, and it may therefore be 
neglected Therefore P,Q is parallel and equal to PP. 

Hence^ the projection of the whole disturbance Pg on any 
straight line, bemg equal to the algebraical sum of the pro- 
jections of PP, and P,(3, will be equal to the algebraical sum 
of the projections of the separate disturbances PP„ PP, 

Now if there are three small distuibing forces^ m, m 
wo may consider the joint action of the two m^, as one 
small disturbing force ; therefore, by what precedes, the total 
disturbance along any axis will be the sum of the separate 
disturbances of m, and of the system m^, but this last is 
the sum of the separate disturbances of and mf, therefore 

B 2 
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the ■whole disturbance equals the sum of the three separate 
disturbances 

This reasoning can evidently be extended to any number 
of foices, and if a:, y, z be the cooidinates of the disturbed 
particle, j>{x^y,z) any function of x, y,®, the disturb<ance pio- 
duced in ^(x,y,s) will be 


where Bx = Sa;, 4- Sajg H- — sum of disturbances along axis of x 

due to sepal ate forces, 

Sy = dy, + Sy, 4 . = • along axis of y, 

Bz = Ss, + 8z^-{- = • along axis of z , 

therefoio B^ (x, y, «) = ^ Sx, + «y, + dz. 


(I(j) 

(lx 


d[/ 




dtf) 

Tz 


= ( r, y, z) + B^<f> (x, y, s) + &c., 


oi total disturbance equals sum of separate disturbances, which 
proves the pioposition. 


4 Since ^(x, y, z) may be the radius vector, or the latitude 
or longitude of the disturbed body, it follows that the total 
disturbance in any of these elements is the sum of the partial 
disturbances. 

Thorofoic m determining the motion of a secondary relative 
to its piimiu'y, as in the picsent case of the moon about the 
earth, where the disturbing effects produced by the sun and 
planets arc small, wc may consider them one at a time, and 
hence the famous pioblein of the Thai Bodies. 

The planets being small and distant, their effect on the 
motion of the moon will not be of sufficient intensity to affect 
the order of approximation to which it is intended to carry 
the solution in the following pages, and our problem is reduced 
to tlio consideration of the three bodies, the sun, earth, and 



ATTRACTIONS. 


5 


moon, acting on one another according to the law of nnivex'sal 
gravitation 

5. But we must still prove another proposition, without 
which the problem would scarcely, though reduced to 
bodies, be less complicated than m its most general form- 

Newton s law refers to partteles^ wheieas the sun, earth, £md 
moon are large spheiical bodies, and it becomes necessary' *0 
examine the mutual action of such bodies Now, it 
that with this law of force, the atti action of one sphero on 
another can be correctly obtained, and leaves the questioxx in 
exactly the same state as if they were particles (Pr'iiitcij) 
lib. I. piop. 75 ) 


Attractions of Spherical Bodies 

6. Let P (fig. 2) be a particle situated at a distance OP == a 
fiom the centie of a xmtform attracting sphere whose density isp 
and ladius OA = c a> c, the particle being without the splioi'c. 

Lot the whole sphere be divided into circular laminae By 
planes pcrpoiidiculai to OP. Let one of these PS == j", 

PQ = and thickness of lamina = Sj? 

Next, let this lamina be divided into concentric iings. 

Jth = r be the ladius of one of these rings and Br its bretvcl'Lli, 

theicforo r = x tan 

B7 = X sqc^O.BO 


The attraction of an element P of this ring on the particle 

^ „ - mass of element , n-n , 

Jr Will be rSxii along jPZi, and the resolved pai't of 


PIP 

this along PO will be 


mass of element 
x^ 8ec^0 


cos 0. 


But the icsultant attraction of the whole iing will clea.x*ly 
bo the sum of the resolved parts along PO of the attractions 
of its constituent elements; therefore, 

attraction of iing = cos 0 — 27rp sm0 ; 

® X sec ^ ^ ’ 
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rcos"i- 

therefore, attraction of whole lamina SQ == 27rpSx I "siii^ dO 

= 'iirphx ^1 - ^ 


Again, = cc*'* -{- — (a — xY = 2ax — (a'* — c ^) ; 

therefore = aSx^ 


and attraction of lamina 


_ fzSz 

= 2,rp(— - 


zSz iS'* + a"* — 


2a'* 




)• 


... ^ I, 1 T o {d^-^C^)z 

attraction of whole sphere = 2Trp ^ — 2 a^ 


|2a Ga'* 
(from z = a — do z=^a-\‘c) 

_ 47rpc* 

~ Ba^ 

M 




where M = mass of sphere = - 

o 

Hence, the attraction of the whole sphere is pjeciscly the same 
as %f the whole 7)ia$& ivere condensed into %ts cenfte 


Cor 1 The attraction of a shell radius = c and thickness 
Sc will he obtained from the preceding cxpicssion by differen- 
tiating it with respect to c, and is 


atti action of shell = ~ — ^ 
a 


mass of shell 


Cor 2 Theicfoie, the atf) action of a heteiogcncous sphere 
on an external particle will he the smne as if the whole mass 
were condensed mto its centre^ provided the dots it y he the same 
at all points equally distant from the cenfie^ for then the whole 
sphere may be considered as the aggregate of an infinite number 
of uniform shells, and by Cor. 1, each acts as if condensed mto 
its centre 
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7. Let us now consider the case of one spheie atti-actinff 
another Suppose P in the preceding article to he an elementary 
particle of a sphere if', whose centre O' suppose at a distance 
a horn 0. Then, since action and reaction are equal and 
opposite, P will attract the whole sphere M just as it would 
do a particle of mass if placed at 0; the same is true of all 
the elementaiy particles which compose the spheie if', therefore 
the sphere if' will attiact the sphere if as if the whole mass 
of the latter were condensed into its centi-e 0 but the attraction 
of the sphere if' on a particle 0 is the same as if the attracting 
sphere were condensed mto its centte O' ; therefore. 

Two spheres attract one another as tf the whole matter of each 
sphere were colhoted at its centre. 

8 This remarkable result, which, as may be shewn, holds 
only when the law of attraction is that of the inverse square 
of the distance, or that of the direct distance, or a combination 
of these by addition or subtraction, reduces the problem of the 
sun, earth, and moon to that of three particles; — the slight 
enor due to the bodies not being perfect spheres will here be 
neglected, being of an order higher than that to winch we 
intend to cany the piesent investigation; this error however 
though very small, is appreciable, and if a nearer approximation 
were required, it would be necessary to have regard to this 
circumstance (See Appendix, Art 100) 
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MOTION RELATIVE TO THE EARTH 

9 When a number of particles aie m motion imder their 
mutual attractions or other foiccs, and the motion relatively to 
one of them is requiied, wc must bring that one to rest, and 
then keep it at lest without alteiing the relative motions of the 
otheis with icspect to it 

Now, fiistly, the chosen particle will be brought to lest by 
giving it at any instant a velocity equal and opposite to that 
which it has at that instant, secondly, it will be kept at rest 
by applying to it acceleiating forces equal and opposite to those 
which act upon it 

Tlicicfoic give the same velocity and apply the same accele- 
rating foices to all the bodies ot the system, and the absolute 
motions about the choHon body, which is now at lest, will be 
the same as their lelativc motions pieviously. 

Prohhm of Two l^oihca 

10, As the sun distuibs tli(5 moon’s motion with respect to 
the eailh, it is impoii:ant to know what that motion would have 
been if no disturbance had ovisted, ov generally • — 

Two hodi(\s att) acting one amtlie^ with fotecs varying directly 
as the ma>ss and in verse fi/ as the stjuaie of the distance^ to deter-- 
mt?w the oihit of one relatively to the other* 

Let ilf, be the masses of the bodies, r the distance between 
them at any time M' being the body whose motion relatively 
to M IS rcqiiued. 

M 

The accelerating force of M on M' equals -p acting towards 
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M' 

while that of M on M equals m the opposite direction 

Theiefore, by the principle above stated, we must apply to both 
M and M’ accelerating forces equal and opposite to this latter 
foice, and M' will mo\e about Affixed, the accelerating force 

on Af' being — — = ijlu\ if yit = Af + Af' and ^ ~ 


Hence, 


d^u 
dO 

wheie }i = r^ ^ ^ twice the area desciibed in a unit of time, 
at 


w = {I4 6cos(0-a)}, 


e and a being constants to be determined by the ciicumstances 
of the motion at any given time 

This IS the equation to a conic section referred to its focus, 

the eccentricity being the semi-latus lectum — , and the angle 

made by the apse line with the prime radius a 

In the relative motion of the moon, or in that of the sun 
about the earth, the orbit would, as observation informs us, be 
an ellipse with small eccentricity, that of the moon being about 
^ and that of the sun ^ 


11. The angle 6 — a between the radius vector and the 
apse line is called the true anomaly 

If n IS the angular velocity of a radius vector which moving 
uniformly would accomplish its i evolution in the same time as 
the true one, both passing through the apse at the same instant, 
then 4 £ — a is called the mean anomaly^ £ being a constant 
depending on the mstant when the body is at the apse, its value 
being also equal to the angle between the prime radius and the 
uniformly revolving one when ^ = 0. 

Thus, if AfT (fig 3) be the fixed line or prime ladius, 

A the apse, 

Af ' the moving body at time i, 

Af/xr the uniformly revolving radius at same time, the 
direction of motion being represented by the arrow 
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Let MD be the position of when t = 0, 
then ^MD = s and is called the epoch,^ 

DMfjb = nt^ 

TMA = a = longitude of the apse ; 
therefore, mean anomaly ~ AMju. — nt + s - 
true auomaly = AMM' = TMM' - TMA - 6 - a 

12. To express the mean anomaly in terms of the true in a 
s&nes ascending according to the powers ofej as far as 

^ _ Stt 2 area 

periodic time ” ”” 

27rh __ h 
27rab ~~ a^ — e^) ’ 
therefore h = nd^ (1 - 

^ ~ 1 
dO h h {1 + 6cos(0- a)}'* 

= ^ (1 - {1 + e cos(0 - a)}"** 

= J (1 - fe') {1 - 2e cos(0- a) + 3e* CQa\e-a)] 

= ^ {1 — 2e cos(^ — a) + |e“ cos2 (^ — a)] ; 

therefore nt->r b = 6 -2e sm(0 - a) + fe'‘ sm2 (0 - a), 
or [nt+e — a] = (^-a) — 2e sin(^ — a) + fe'* sm2(^ — a), 
the required relation. 

13. To eatress the true anomaly in terms of the mean to the 
same order of a]y^oxvmciA%on 

5 - « = M# + e - a + 2e sin (^ - a) - le* sin 2 (^ - a) . (1) ; 

.. d-a — nt + e — a first approximation 

^ The mtroductioxi of the epoch is avoided in the Lunar theorj by a par- 
ticular assumption (Art 34) , but m the Planetary it forms one ot the elements 
of the 01 bit 
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bubstituting this m the first small terms of (1), we get 
^ — a = + £ — a -h 2e sm — a) a second approximation. 

Substitute the second approximation in that small term of (1) 
which IS multiplied by e, and the fiist approximation in that 
multiplied by e*'*, the result will be correct to that tenuj and 
gives 

0 — a = nt~j~e-^a’i‘2e sinfwi + e — a + 2e sin(n^ + s — a)} 

- fe'* sm2te + £-a) 

= nt + e — a -f 2e sin (wif + £ — a) 

+ 4<3‘^cos(w^ 4- £ — a) sm(^^^ 4 e — a) — ■|e‘* sin 2 (ni + e — a) 

= nt -{- e — a 4 26 sin (wif 4 £ — a) 4- sin2 4 e — a), 

the required relation. 

The development could be carried on by the same process 
to any power of 6, but in what follows we shall not require 
it further than 


Problem of Three Bodies 

14. In order to fix the position of the moon with respect 
to the centre of the earth, which, by means of the process 
described m Art (9), is reduced to and kept at rest, we must 
have some determinate invariable plane passing through the 
earth’s centre to which the motion may be referred 

The plane which passes through the earth’s centie and the 
direction of the sun’s motion at any instant is called the true 
ecliptic, and as a first rough result of calculation, obtamed on 
supposition of the sun and earth bemg the only bodies m the 
universe, this plane, m which, according to the last section, an 
elliptic orbit would be described, is a fixed plane but this is no 
longer the case when we take into account the disturbances 
produced by the moon and planets, and it becomes necessary 
to substitute some other plane of reference unaffected by these 
disturbances 
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Theory teaches us that such a plane exists,* but as its deter- 
mination can only be the work of time, the following theorem 
will supply us with a plane whose motion is extremely slow, 
and it may for a long penod and to a degiee of approximation 
far beyond that to which we shall carry our investigations, be 
considered as fixed and coincidmg with its position at 
present 

15 The centre of gravity of the earth and moon describes 
relatively to the sun an orbit very nearly in one plane and 
elliptic y the square of the ratio of the distances of the moon 
and sun from the earth hemg neglected f 

E 


Let /?, M be the centres of the sun, eailh, and moon, 
O the centre of gravity of the two last Now” the motion of 
G IS the same as if the whole mass E’\-M were collected there 
and acted on by forces equal and paiallel to the moving foices 
which act on E and M The whole force on G is therefoie 
m the plane 8EM^ join 80 



* See Pomsotf “ Theorie et determination de Veqxmteur du sysUme solaire,** 
where lie proves that an invariable plane exists for the solar system, that is, 
a plane whose position relatively to the fixed stars will always be the same 
whatever changes the orbits of the planets may experience , but as its position 
depends on the moments of mertia of the sun, planets, and satelhtes, and 
therefore on their internal conformation, it cannot be determined a piion^ 
and ages must elapse before observation can furnish sufiioient data for domg 
so d posteriori 

This result Poinsot obtains on the supposition that the solar system is 
a free system , but it is possible, as he furthermore remarks, nay probable, 
that the stars exert some action upon it, it follows that this invariable piano 
may itself be variable, though the change must, accordmg to our ideas of 
time and space, be indefinitely slow and small. 

t This ratio is about tov and, as we shall see Art (21), such a quantity 
we shall consider as of the 2nd order of small quantities, and its square there- 
fore of the 4th order Our investigations are earned to the 2nd order only 
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Let L8Cr]!i!L = ct), and let m’ be the stin’s absolute force 

Moving force on jB* = in E8. 

blL 

moving force on If = in MS 


These^ applied to G parallel to themselves are equivalent to 



m'JSGE m’MQM , , 

8E^ /SM " direction GM^ 


and 

m'ESG , m'MSG 

8IP SM" 

G8, 


But 

8E^ = 8a'‘ +GE^ + 2 8G GE cos 

= 8G^ + GM‘ -2 8G GM cos co 


Hence 

1 1 S GE 1 

SE^ ~ 8G^ 8G* 1 

1 1 SG3I 1 

8M^ ~ SG^ SG* i 

1 omitting 

(EM^^ 
[SGJ • 

Now E GE 

= M.GM= {M-^-E) Therefore the ac- 


celeiatmg force in the diiection GM 


3m' GM GE 

— -• COSG) 

= — 3 (accelerating force of sun on O) cos to 

= 0 accoiding to the standard of approximation adopted. 

Hence the only force on is a cential force tending to 8^ 
and therefore the motion of G will be in one plane 


♦ In strictness it would be necessary, since we haye brought »S to rest, to 
apply to both M and and therefore to G, accelerating forces equal and 
opposite to those which E and M themselves exert on S , but the mass of S is 
so large compared with those of E and Jkf, that we may safely neglect these 
forces in this approximate determination of the path of G, the error being of 


a still higher order than that introduced by the neglect of 




14 


LUNAR THEORY 


Again, the acceleiatmg force in G8 

_m'SG Z{EGE^MQ:M) ) 

m+e\ SG^ ■ 

= to the same appioximation, 

Theiefore the orbit of Q about 8 is very approximately an 
ellipse with 8 in the focus, and the plane of this ellipse is, as 
far as our mvestigations are concerned, a fixed plane when 
8 IS fixed 

This fixed plane is called the plane of the ecliptic^ or simply 
the ecliptic 

16 A plane through the eaith’s centre paiallel to the 
ecliptic will be the plane of reference we require (14) and will 
become a fixed plane when we bring the earth’s centre to rest, 
the ecliptic then making small monthly oscillations from one 
side to the other of our fixed plane 

17. The sun will have a latitude always of the same name 
as that of the moon, and deducible from it, when E8^ EM^ 
and the ratio of the masses of the earth and moon are known. 
For if 8EM' be this fixed plane through E (fig. 4), /S', G\ M\ 
the projections of /S', (?, if, 

then, sin(sun’s lat ) 

— QTTQ^ __ sin(moon’slat ) 

M EM . , , . . 

;^sm(moon’slat) 

Now, from observation it is known that M is about of E^ 

and EM EjS^ 

therefore sm (sun’s lat.) = ^ ) nearly 

‘ 32400 

And as the moon’s latitude never exceeds 5° 9', the sun’s lati- 
tude will always be less than 1" 



ERRORS IN sun’s POSITION 


15 


Again, with respect to the sun’s longitude let ET be the 
direction of the fiist point of Aries, — that is, a fixed line in 
our plane of reference fiom which the longitudes of the bodies 
are reckoned = ff the sun’s longitude 

The diffeience in the sun’s longitude, as seen from E or from 
will be the angle ES'G’, 

EG' 

s,mE8'G'=^ sm© = approximately, if ',8' = 

theiefore sxaES' Q' nevei exceeds 
tlierefore E8'G' is a small angle not exceeding 7" 

Also ES' - 8' G'< EG' < 8’ G'. 

Now, by assuming the longitude and distance of tbe sun 
as seen from A? to be tbe same as when seen from (7, we 
commit the above small errois m the position of 8, that is, 
we assume the sun to be at 8" mstead of N, 8' 8" being drawn 
equal and parallel to G'E If our object were the deter- 
mination of the sun’s position, it would be necessary to take 
this into account, but the consequent small errors introduced 
in the disturbance of the moon will cleaily, on account of 
the great distance of the sun, be of a far higher order, and 
may therefore be neglected. 

18. Hence we may assume that the motion of the sun 
about the earth at rest is an ellipse havmg the earth for its 
focus, and its equation 

u' = a' (1 -l-e' co8(0'— f)}, 

and we aio safe tliat no appreciable error will ensue in tlio 
determination of the moon’s place* 


* 7 hat IS, as far as the tlirao bodies alono ore concerned but, since the 
attractions of the planets may, and in fact do, disturb the elliptic orbit of the 
sun about G, tlio same cause will disturb the aimmcd oibit about E. A re- 
markable result of this disturbance is noticed m Appendix, Art (99) 
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CHAPTEK III 


BIGOKOUS DIPPERENTIAIi EQUATIONS OP THE MOON’S MOTION 
AND APPROXIMATE EXPRESSIONS OP THE POECES 


19 The eaith having been i educed to rest by the piocess 
described m Ait (9), its centre may be taken as origin of 
coordinates, the fixed plane of icfercncc as plane of xy and the 
hue A'Y as axis of x (fig 5) 

Lot 9 , 5 be the cooidinatcs of the pioiection ilf' of the moon 
on the plane xy, s the tanyent of the moon’s latitude MJEW. 
Also let the acceloiating forces which act on the moon be 
resolved into these tlnoo . 

P paralkd to the projected ladius ilf and towaxh the caith, 

T paiallol to the plane xy, pci-pendicular to P and in the direc- 
tion of 9 tuermsmg, 

8 perpendicular to the fixed plane and towards it 

d^x p® „,?/ 


whence 


also 


dt r r 


(Py 

' dt* 






Tr, 


d‘x , d*y p 

+y-j^i 


dt* 




= - 8 . 
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And introducing polar cooidmates, these equations become 


— = 
dt V dt) 


Tr 


d^ \dt) 


d^irs) 

~W~ 


= - B.... 


• • • W) 

(ll), 

. .. (m) 


20 These three equations for determining the moon’s 
motion take the time t for independent variable, but it will 
be more convenient m the following process to consider the 
longitude as such, and oui next step will be to change the 
mdependent variable from t to 6. 

From (i) we get 



d_ / a _ 27^,^. 

dt dt\ dt) ~ dt^ 


therefoie 

{t‘ + 2| IVW0, 



(in 

= suppose, whence II 

= 

Thcrefoie 

dB K „ i f 1 



dt \ 1 



dd Hu‘ , „ //, ^ r t 


Again, 

dr dr dd du 

dt~~^ dd' dt ^ dd^ 

d^r d ( du\ ^ a d / dri\ 


Substitute 

(“)» 


therefore 



or 



■(A); 


0 
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aad writing for IP and for ^ their values + 2 and ^ 
the equation becomes ^ ^ 


<Pu 


P 

AV 


T 


dd 


1 + 2 


r T 

}h‘u' 


ie 


(/?} 


This 18 ^lled the differential equation of the moon's radius vectoi 
Liastly, 


- 8 = 


from (A), 
therefore 

therefore 



— \ rr ( „ ds du 


dt 


dd 

d^u 


dii 

Te. 

_ du 

''de. 


Ps = EM 


5-^, + MS 


I IT 2 




■ps 

de^ 


+ S = 


fd^s 

Ps- 8 
Av ■ 


dd 

dd ' 

_T_ ^ 

AV dd 


1+2 


Jav 


:dd 


ij) 


This IS the differential equation of the moon's latitude 
If the three equations (a), (^), (^) could be integrated undei 
these general forms, then, since they are perfectly rigorous, the 
problem of the moon’s motion would be completely solved for 
as 0 y fom vaiiaUes m, d, s, and t are mvolved (the acceleraW 
Ws P, T, and 8 are functions of these), the values of three 
ot them, as m, d, s, could be obtamed in teims of the fourth f 
that is, the radius vector, longitude, and latitude would be 
known correspondmg to a given tune. 

r Ilie mtegration has never yet been effected, except 

for particular values of P, T, and 8; and the method which w^s 
are m ^nsequence forced to adopt, is tliat of successive ap- 
proximation, by which the values of .m, d, and s are obtamed 
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ORDERS OF j^iMALL QUANTITIES 

in a series, the terms pioccedmg accoidmg to ascendmg poweis 
of small fi actions, some one being chosen as a standaid with 
which all othcis aic compared, and the order of the appiox- 
imation is estimated by the highest power of the small fi actions 
retained. 

It IS usual to consider as a small fraction of the first oidei, 

consequently of ^ is second 

is thud 

and so on, othei fractions being considered as of the 1st, 2nd, &c 
orders, accoidmg as they more neaily coincide with &c 

22 It is ncccssaiy theiefore, before we can approximate 
at all, that we should have a previous knowledge (a rough one 
is sufhcient) of the values of some of the quantities mvolved in 
oui mvostigations ; and foi this knowledge we must have le- 
couiso to ohscivation. 

We shall thcicfoie assume as data tlie following results of 
ohscivation ; 

(1) The moon moves in longitude about thirteen tunes as 

fast as the sun The accuiate ratio of the mea 7 i motions in 
longitude repicscnted by m is theiefore about be 

consideicd as of the 1st order * 

(2) The sun’s distance from the earth is about 400 tunes 
as great as the moon’s distance. 

Hcncc the latio of the mean distances = is of the second 
order- f 

(3) The cccentncity e' of the elliptic orbit which the sun 
appioximatoly describes about the earth is about and this, 
appi caching nearer m value to than to will be con- 
flidorod as of the 1st order 


* TIub apiiroxxmate yaluo of m is easily obtained , — the moon is found to 
perfonn tho tour of the heavens, returning to the some position among the 
fixod stars, m about 271 4tiys , the sun takes 365i days to accomphsh the same 
jouniey 

t Tho distances of the luminaries may be calculated from their horizontal 
parallaxes, found by obscivations made at remote geographical stations 

C2 
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( ) uiing one revolution, the moon moves pretty accurately 
map ane mclined to the plane of the ecliptic at an angle whose 
tangent is about and therefore of the 1st order * 

(5) Its orbit in this plane is very nearly an ellipse having 
e centre of the earth in its focus, and whose eccentricity is 
about equal to our standard of small fractions of the 1st order, 
VIZ. ^ ; and this will also he very nearly true of the proiection 
ot the orhit on the plane of the ecliptic-f 


To calculate the values of P, T, 8 

23 We are now in possession of the data requisite foi 
beginning our approximations, and we shall proceed to the 
determination of the values of P, T, and S m terms of the co- 
ordinates of the positions of the sun and moon. 

Let P, E, M (fig. 6) be the centres of the sun, earth, and 
moon, ’ ’ 

m', M their masses, 

-E', if', the projections on the plane of the ecliptic, 

O the centre of gravity of E and if, 

fi — E + M, 

MGM = tan 's = moon’s latitude, 

E'T the direction of the first point of Aries, 

8G = r =— , L HE' S= P = longitude of sun. 


distances being then dravn firom a ^ T P^Portional to the 

trenubesmaih out the form of the moon's orbit direotionB, the ex- 
sun wlT detennmation of the eccentricity of ihe 


exphessionjs of the foeoeis 
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M'E* — r r= 1 j ^ TE'M* = 0 = longitude of moon, 

-• SE'M' = 0 — 0' = difference of longitude of sun and moon. 

The forces we have to take mto account are, according to 
Alt (9), the forces which act directly on and forces equal 
and opposite to those which act on E ^ — ^these last being apphed 
to the whole system so that E may be a fixed pomt 

Atti action of S upon M = m -MiS, eqm valent to 


m 


m 


, MG 
SM^' 
t ' 


attraction of E upon M = 


E 

ME^ 


m MG^ 

parallel to G8\ 
mME^ 


atti action of 8 upon E = m E8^ equivalent to 

{ ,EG 


in EG, 


[m' • parallel to G8'j 

M 

attraction of M upon E = in EM. 


yiicrefore, the whole attraction upon M^ when E is brought 
to ICSt, IS 

E^M ,fMQEG\ „„ 

MB^ + ^ [sip ^ SB) “ 

and mV parallel to GS. 


These expressions of the accelerating forces on M are rigorous, 
and can be expressed m terms of the masses and coordinates of 
the bodies ; but since our mvestigations will be carried only to 
the second older, it will be sufficient if, in the preceding, we 
neglect small quantities of the fourth and higher oideis. 
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Now, 8M^ = 8M''‘ + MU'* 

= 8G^ + GM'^ - 28G GM‘ cos80M‘ + MM'^ 


therefore 


= r'^ ("i _ 2 C0S8GM' + ^ 
\ r 

1 „ 1 f, , %GM' 

8M^~ r'^ 1^ ~ 


for 6-6' ox 8E'M' differs from 8GM' by less than 7", Art (17), 

/ Q-M\^ 

aiid ( IS neglected, being of the fourth older, Ait. (22) 


1 1 


Similarly, ^ ^ cos {6 - 0')| i 

therefore, the acceleiatmg forces on the moon are approximately 
{MG + GE) in direction ME, 

and ( GM' + GE') co&{6 — 6') parallel to ; 

whence P = cosMGM - ^ ME' cos'‘(0- 6') 

At mV 3mV . 

= T77Ti'^^~'^ 

r (1 + sy 

= /im“ (1-1 s‘) - {^ + f cos2(0 - 61')}, 

fKfiW 

P = - ^ M'E oob{ 6-6') sm{d- 6') 

= - f sm2(0 - 6') 

^ = [-Me‘ + Y^^^)^mgm' 

_ f /* to'? V(l+g^J ) 3 

V(l + a")'^ »*'•’ 1^(1 + sO 


r'^ jV(l + sO 


2 / ft *1\ I ^ ® 

= /iM(s-|-s ) + -—-; 

w 


Ps- 8 = 


+ ^cos2(0-0') 



OKDER OF DISTURBING FORCE. 
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24 The differential equations in Art (20), when these 
values of the forces are substituted m them, would contain a 
new variable 6\ but we shall find means to establish a con- 
nexion between 6^ and ff which will enable us to eliminate 6\ 
They will, however, be still incapable of solution except 
by successive approximation j but before proceeding to this, it 
will be important to consider the order of the disturhing effect 
of the sun’s action, compared with the direct action of the earth 
Now, if we examme the values of P, T, and 8^ it will be found 
that the most important of the terms containing m', which are 
clearly the disturbmg forces smce they depend upon the sun, 

are involved in the form ^ , while those mdependent of the 
sun’s action enter m the form ~ . 

We must therefore find the order of ^ compaied with ^ , 


or of 


m 


^3 


r 7" 

Now the orbits being nearly circular, and m the ratio of the 
mean motions. Art (22), we have 

mean motion of sun _ periodic time of moon 
mean motion of moon “ peiiodic time of sun 


7n = 


27rr* 27rr' 


A 


m 


theiefore 


i ’ 

7n' 

r' 




or the disturbmg force of the sun is of the second order. 


t 
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CHAPTER IV 


INTEGEATION OP THE DIPPEEENTIAL EQUATION S. 


Section I 


Q&nefdl process desc'nhed. 


25. The differential equations wkeh ve have obtained are, 
as already stated, incapable of solution m then- general forms- 
and even when P, T, &c have been replaced by their values’, 
the mtegration cannot be effected, and we must proceed by 
successive approximation, 

Iirstljj neglect the disturbing force of the sun which is of 
the second order, and also the moon’s latitude, which, as will 
be seen by referrmg to the expressions for the forces (23), will 
either enter to the second power or else m combmation with 
the disturbing force. 

When this IS done the equations can be mtegrated, and 
values of M and s obtained in terms of 6 correct to the same 
order of approximation as the differential equations themselves, 
that IS, to the first order, and this value of u will enable us 
to get the connexion between 0 and < to the same order 


hoyeyer, bear m ramd that the equations thus mtegrated are 
equations of the moon’s motions, but only approximate 

nThel. ’ “ “ that the results obSned may 

not be even approximate forms of the tme soluUons 

f ’’y compaiing them with 

^hat we aheady knov of the motion from obscnaUon. Ld this picvious 
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knowledge, in the event of then not being approximations, will probably 
suggest such modifications of them as will lender them so *] 

The integration of the equations (a), (/3), ( 7 ) can be pei- 
formed when the second members are circular functions of ff] 
and as the fiist approximation will give us the values of u and s 
m that form, when these values are admissible and carried into 
the expressions for the foices, they also will be expressed as 
functions of and we can proceed to a higher approximation. 

The ,new approximate values of P, T, 8 are then made use 
of to reduce the second members of the differential equations 
to functions of retainmg those terms of the expressions which 
are of the second order 

The equations are agam mtegrable, and this being done, 
the values of w, 5 , t will be obtained correctly to the second 
order These values introduced m the same manner in the 
second members, and terms of the next higher order retained, 
will lead to a third approximation, and so on, to any order, 
except that if we wish to carry it on beyond the third, the 
approximate values of the forces, given in Art (23), would no 
longer he sufficiently exact 

26. There is, however, a peculiarity in these equations, when 
solved by this process, which we must notice We have said 
that to obtain the values to any order, all terms up to that 
order must be retained in the second members but it may 
happen that a term of an order beyond that to which we are 
working would, if retained, be so altered by the integration 
as to come withm the proposed order. 

Such teims must therefore not be rejected, and we shall 
proceed to examme by what means they may be recognised. 


* We cannot say a priott tliat tlie solution of n = f {u, 6) shall be of 
(Pzt 

tlie same foim as that of oven though should be a very 

appioximate value of / (?«, 0), but there is a piesumption in favoui of such a 
supjpositioii 
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27, Suppose then that after an approximation to a ceitain 
Older, the substitutions for the next step have brought the 
equation m w to the foim 

+ ^COs(^0+^)+ .... , 

where the coefficient O is one order beyond that which we 
intend to retain The solution of this equation will be of the 
form 

u— cos(_p0 + ^) -f 

G* bemg a constant to be detcrmmed by putting this value of 
u m the differential equation, 

G 

■whence G' = j , 

1 -p-*’ 

from which we learn that if p differs veiy little from 1, G' will 
ho one order lower than G^ and will come within om* proposed 
appioximation, and consequently the term O cos(j?0 + il) must 
be retained in the differential equation 

The equation of the moon’s latitude being of the same form 
as that of the radius vectoi, the same remarks apply to it 


28. Agam, m finding the connexion between the longitude 
and the time (one of the principal objects of the Theory), we 
must use equation (a), Art. (20), 

^ 1 


Now, havmg developed the second member and substituted 
foi M, &c. their values m terms of let it become 

'^0 “ ”h ^ cos^^^ jKj 4“ , , , » , y 

hence, t=s + ^ sm{q6 + B) + . , . 

Therefore, when a is of the first order, — will be one order lower 

2 

than and the toim will have risen in impoitancc by the 
integiation 
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But yet further, if such teims occur in -3 , they will be twice 

f T 

increased m value; for they increase once ui forming dQ^ 
and once again, as above, in finding U 

T 

Since such terms occur m the development of and also 

of ~ , on account of their previously bemg found in w, we must 

examine how they appear in the ditferential equation that gives 
that we may recognise and ret am them at the outset. Now, 
by referring to the last article, we see that when p is very small 

1 T 

G and G will be of the same order ; and in -1 , -o the order 

of the term will still be the same ; so that all the teims which 
it will be necessary to retain are known at the outset. 


29. We have, therefore, the folio wmg rule : 

In approximating to any given order ^ we mnstj %n the diffe-- 
rential equations for u and s, retain periodical terms ONE ORDER 
beyond the proposed one^ when the coefficient of 6 in their argument 
is nearly equal to 1 or 0^ and terms in which the coefficient of 
the argument is nearly equal to 0, must be retained TWO ORDERS 


higlm* than the proposed approximation when they occur in . 

If wo wished to obtain u only and not f, there would be no 
necessity for letammg those terms of a more advanced older 
111 which the coefficient of 6 nearly equals 0. 


Section II. 

To solve the Equations to the first order 

30. We shall m this first step neglect the terms which 
depend on the disturbmg force, i e those terms which contain 
wi', for we have seen, Ait. (24), that such terms will be of the 
second 01 del 
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The differential equations may be written under the moie 
convenient form 


d^u _ P 

d^s ^ _Fs-S 

JiW 


T du /d^u 
d^s 


AV dd 
T & 
AV dd 




dff‘ 


+ 




The latitude s of the moon can never exceed the inclination 
of the orbit to the ecliptic , but this mclmation is of the first 
Older, therefore s is at least of the first order and s‘ may be 
neglected 

Therefore, from Aid ( 23 ), 





and tho differential equations become 


d^u fju 


d^s 


5 = 0 ; 


whence ^ = p {1 + c cos(0 — a)}, or, writing a for ~ , 

u = a {1 4 e cos(i 9 -a)} (?7J, 

and 5 = A sm ~ 7) . ( ;SJ, 

C5 a, yfc, 7 being the foui constants introduced by integration 


31 , These results are in perfect agreement with what ^ough 
observations had already taught us concemmg the moon’s motion 
Art, (22) , for 

w = a {1 4 5 cos(0 — a)} 

represents motion in an ellipse about the earth as focus. 

Again, 5 = & sin (0 — 7) indicates motion m a plane inclined 
to the ecliptic at an angle tan“^A. 

For, if TOifef' be the ecliptic, (fig 7) 

M the moon’s place, 

MM' an aic perpendicular to tho ecliptic, 

Tiir = <9, 


then 
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and if TO be taken equal to y, and OM^ joined by an are of 
groat circle, we have 

sin OM! = tan JO/' cotMOM' ; 

01 sin (5 — 7) = s cotMOM', 

which, compared with the equation above, shews that 
MOM' = tan'^Jj 

Therefoie, the moon is in a plane passing through a hxed point 
0 and making a constant angle with the ecliptic ; or, the moon 
moves in a plane. 

32 What the equations can not teach us, however, and for 
which we must have recouise to our observations, is the approx- 
imate magnitude of the quantities e and h By referring to 
Art. (22), we see that e is about ^ and k about -j^, that is, 
both these quantities are of the first order Their exact values 
cannot yet bo obtained the means of domg so from multiplied 
observations will be mdicated further on. 

The values of a and 7 introduced in the above solutions are 
respectively the longitude of tlie apse and of the node 

33 . Lastly, to find the connexion between t and 6 , the 
equation (a) becomes, makmg T= 0 , 

^ __ 1 _ _ _1 1 

dd hu^ ha‘ {1 + e cos(0 — «)}■' 

Now this IS the very same equation that we had connecting 
t and 6 in the problem of two lodws, Art. (12), as we ought to 
expect, since we have neglected the sun’s action Therefore, 
if p be the moon’s mean angular velocity, we should, following 
the same process as in the article referred to, arrive at the 
result 

6 + e + 2e sm(^« + e — a) + sm2 {pt + B — d) + , 

which IS correct only to the first order, since we have rejected 
some terms of the second order by neglectmg the disturbing 
force 
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34 The aibitiaiy constant e, introduced in the pioccss of 
integration, can be got rid of by a proper assumption this 
assumption is, that tire time t is reckoned fiom the instant when 
the mean value of 0 is zeio * 

Foi, smce the mean value of 0, fomid by lejeetmg the 
peiiodical teims, is^jt + s; if, when this vanishes, < = 0, we 
must have s = 0; theicfoie 

0 = pt + 'Ze sin [p — a) 0^, 

correct to the first ordei f 


35 We have now obtained three lesults, Z7j, S',, @„ 
as solutions to the fiist older of oui diffciential equations, 
and we must employ them to obtain the next appioximatc 
solutions . but before 27, and 8^ can bo so employed they must 
be slightly modified, in such a maimci liowevei as not to intei- 
fere with then degree of appioximation 

The necessity foi such a modification will appear from tlio 
following considerations . 

' Suppose we proceed witli the values alicady obtained; wc 
have, by Ait. (23), 


P 


a 

¥ 


(1- K) - 


mu 


&c . 




mu 

Wd 


{1 -he cos (0- a)}"*, 


+ A. cos (0 — a) -h . j 


* ‘When a function of a variable contams periodical ternis winch go through 
all their changes positive and negative as the variable increases continuously, 
the mean mlm of the function is the part which is independent of the periodical 
tenns 

t We shall also employ this method of correcting the integral in our next 
approximation to the value of 6 in terms of t , and if wo purposed to carry oui 
approximations to a higher order than the second, we should still adopt the 
same value, that is, zero, for tlie arbitr^uy constant mtioduced by the integra- 
tion To shew the advantage of thus correctmg with icspcct to mean values 
suppose we lechoned the time from some dejvnite value of 0, for instance when 
08=0 , then, in the first appioximation, 

O=e + 20 sm(a — a). 
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or substituting for from Art (45), 

f 1 — f cos(c0— a) — cos [( 2 — 2 m) 0 — 2 ^} 

^+M=a, + ^m!‘e cos {{2 — 2 m-c )0 — 2 ^ + a.} 

diw 

-f cos 2 (^0—7) — f wV cos(m0+ jS — f); 

therefore, 

’1 — 1^2 — ^992-^ + Q cos(c0-“ a) + cos{(2 — 2m) 0 — 2^} 
=a, + cos{(2 — 27/^ — c) 0 — 2/3 + a} 

— \K^Q.o^^[gd—y)- §mVcos(m0+/3- 5) 
If we compare this with the value of u found Art (48), we see 
that the elliptic inequality, the reduction, and the annual equation 
are due to the central or radial force, as also one half of the 
variation and about a thud of the evection 

It would perhaps be proper to separate the absolute cential 
foice from the cential disturling force, the terms due to the 
latter are those which contam m, theiefore, the elliptic in- 
equality and the reduction are the effects of the formei, except 
that m the elliptic inequality the introduction of c, or the motion 
of the apse, is due to the disturbing force 

(2) To detmm%ne the effect of the tangential disturbing force 
Let the cential distmbzng force be zero , 

then 77-7 = (1 - f5^) = a neglecting the inclmation, 

tl K/ fl 

^8 = - sm{(2 - 2 m) 9 - 2/S} 

+ 3?n^e sin[(2 — 2m — c) 0 — 2yS + a] 
omittmg the term of the fourth order, 

u — a{l + e cos(c0 — a)} ; 

therefore ~ cos{(2 — 2m — c) 0 — 2;S + a}. 


tZ0 == cos{(2 — 2m) 0 - 2y8} 

- 3m^e cos{(2 — 2m - c) 0 - 2/3 -h a} , 


— + 9/ = a to the fiist oidei 



76 


LUNAR TIIEORT! 


Substituting these values in 

d^u P T ^ 

hV dd 


the differential equation 



whence 


= a 


u = a 


1 _ cos{(2 - 2m) 0 - 2/3} 

+ ^m'^e cos{(2 - 2m — o) 0 — 2/3 + a] , 

1 + e cos(o0 — a) + cor [(2 - 2///) 0 — 2^} 
+ cos{(2 - 2wi — c) 0 — 2jS -f a) 


We have here the remaining half of the vaiiatioii and lathoi 
more than two-thirds of the cvcction as tlic clleclh of the t<in- 
gential disturbance Also c=l, oij to the hocoikI oideij the tan- 
gential force has no effect on the motion of the 

The mequahties in the longitude could he easily ohtaiiicd 
fiom the relation 


— — 


I 




1 ’ 


but they would load to the very same eouclusions as the discus- 
sion of the values of u 


To calculate the value of c to the thud otdei 

94 We must hero nuke use of llic icsnlts which the ap- 
proximations to the second oidei Iiavc liunislMsl; but as the 
value of c is determined by that term of the differential equa- 
tion whose aigumcnt is i'0 — a, we need only considei those 
terms which by their combinations will huid to that om* without 
rising to a higher order tlian the hmith. 

We shall simplify the aigumcnts by omitting 0^ which 
can easily be supplied by remarking that - a and in0 — (i 
always enter as one symbol, c and m will theudbie he Mutficient 
to distmgmsh them This only applies to the arguments. 

We have, Aits (48), (23), 

u — a[l e cos(r) 4 //r cos(2 -■ 2///) f C‘Oh(2 — 2;^/ — < )}, 
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P a* 

= a - ^ {1 + 3 cos{2 - 2m)}, 


_|m^“ sm(2-2m) 

Fiom these, we obtain 
P 

= a - {1 H- 3 cos(2 - 2m)} {1 — Be cos(c) 

-f — cos(2 — 2m — c)} 
= a 4 cos(c) + ^m“ae cos(c), 

P 

p -3 = — Im*'* sin(2 - 2m) {1 — 4^ cos(c) — 4m'*cos(2 — 2m) 

It u 

— ^me cos(2 — 2m — c)} 
= — sin(2 - 2m) + 3m®e sm(2 — 2m — c) + ^m®e sm(c), 

T du 

^ = {— f sin(2 — 2m)} { - ^mae 8m(2 — 2m — e)} 

+ {3m^e sm(2 — 2m - c)} {— 2mfa sm(2 - 2m)} 
= f^m^ae cos(o), 

the other term is of the fifth order, 

cos(2 — 2m) — 3m*e cos(2 - 2m - c) — <^m’e cos(c), 

2 + m) = 2a{l + - 3m“ cos(2 - 2m)+ cos(2 - 2m- c) }; 

therefore 2 4- uj dd ^ . — ^m^ae cos(c). 


Substituting in the equation for w, 


d^u 


T du 


(d^u ^ \ fT 


dd 


= a{l + (fm®e+ + ^m®e) cos(c) + .. } 

= a{l + (fm“e+ ajf^m’e) cos(c) + . }. 

Assume u = a}! + e cos(o) + . } ; 
therefore ae(l — c”) = (f m'^e + a ; 

therefore c = 1 — fm'* — ^/m . 



78 


LUNAR THEORY 


To find the vtdm of g to the third order 
95 This IS to be obtained m a very similai manner from the 
equation We shall, in the aigument, wiite g foi 

s = 7tf{sin(9') + |w sin(2 - 2 ot - 9')}, 


Fs-B 


2u^ 


{1 + cos(2 — 2m)}, 


T 3mV 


sm(2 — 2m) 


Whence 

~3 - = - fm% {1 + cos(2 - 2m)} {siii(.<7) + |m sm(2 - 2m - //)} 

fl Vi 

= - lorei (1 - xV«) ™(.7)5 

^ = Tc cos(^) 4 %mh cos(2 - 2m - g) ; 


therefoie 


T ds 

Wdd 

d& 


13 4 s = terms 


- sm},^), 

of the third 01 dor; 


therefore 2 




4 



dd — 0 , to the fourth order. 


Now, 


, Po - ^ T ds 
dB^'^^~ /tV 7iV dB 



= - ’mVc (I -iV»- f\W') 
Assume s == 7». sin(^^) ; 
therefore h{l —g^) = — ]m7». + 


96. Hence, to the third order of approximation, 

mean motion of apse _ 1 “ c. _ f _ 8 + 75m 

mean motion of node 1 8 — 3m ’ 



PARALLACTIC INEQUALITY 


79 


and since m = ^ nearly, we see that the moon’s apse piogredes 
neauy twice as fast as the node regiedes 

In the case of one of Jupiter’s satellites, mis extiemely small, 
for the peiiodic time lound Jupiter is only a few of oui' days, 
and the pcuodic tune of Jupiter round the sun is 12 of oui 
ycais, and therefoie m, the ratio of these periods, is very small. 

Hence, the apse of one of Jupitei’s satelhtes progredes along 
Jupiter s ecliptic, with pretty neaily the same velocity as the 
node legredes, assuming these motions to be due to the sun’s 
distuibmg force, they are, howevei, pimcipally due to the 
oblateness of the planet. 


JParallacUc Inequality 


97, In carrying on the approximations to a higher ordei, 
it IS found, as we stated Art (55), that the expressions for 

jp fjp — - A/* ' 

AV contain the terms - ^ cos(g-g'), and 

E — M a! , n /in -1 

— — sin 0 ) respectively 


Smee ~ = nearly, is of the second ordei, these terms are 

of the fourth order, but the coefficient of 6 being near unity, 
they will become important m w, and therefore in 5, Art (27). 

Wc can easily obtain the terms to which they give rise in 
the values of u and 0, 




= a 


1 + 


E+M 


^ cos{(l-m) , 


T 


E M a 

u = a[l — e cos(c6 — a) 4 . . }, 




^ = ae sm(c0 — a) . 
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therefore 

^ = 0, to the fouith order, 

AV dd ’ 

^ + w = ^5 to the second order, 
d(7 

Snbstitating ia the differentia,! equation for u, we get 
g + « = a[l+.. -V»^^“§^|coa{(l-m)0-^}4 ] 

Assume ii = <2rn- - 


therefore 


therefore u 


1 — [l-m] 




= a 1 + . — M' 


E-Md 
m ^ — Tj.— cos 
E+Jf a 




98. The corresponding term m the value of d wiU also be 
of the third order, 

f T 

f TT-a d0 IS of the fourth order and will not nse m t : 

JAV ’ 

theiefore 

If,. , IS a' ,u NSI oi 1 

* = |[^+ •-)-*^’"f^^sin{(l-^)0-y3] + ..], 

and 6= p +. . . - ^ 8m{(l-w)2J{-y3} +. . 
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This term, whose argument is the angiilai distance of the 
sim and inooiij is called iha pm allactic inequality on account of 
its use in the dctciinuiation of the sun’s parallax^ to which 
puipose it was fiist applied by Mayci by compaiing the ana- 
lytical cxpicssion of this coefficient with its value as deduced 

from obseivation The values of on and of and tliereforcj of 

E-M. , a' 

pi<itty accuiatcly known, ~ will be determined, 

that IS, the ratio of the sun’s paialLix to that of the moon but 
the moon’s parallax is well known, tlieictoie also, that of the 
sun can be calculated The value so obtained for the sun’s 
parallax is 8 63221 ", while those given by the two last transits 
of Venus fall between 8*5" and 8 1”^- 


Soculm Acceki atwn 

99 Halley, by the compai ison of ancient and modern eclipses, 
found that the moon’s mean revolution is now performed in a 
shelter' time than at the epoch of the rccoidcd Chaldean and 
Babylonian eclipses. The explanation of this phenomenon, 
called the spcnhu (ucelpiation of the onoon\ mean onoUon^ was 
for a long time unknown it was at last satisfactorily given by 
Laplace. 

The value of Art ( 50 ), on which the length of the mean 
peiiod depends, is found, when the appi oxiination is caiiicd to 
a higher cider, to contain tlic quantity e the eccentricity of tlio 
eai til’s orbit Now, this ccccntiicity m undergoing a slow but 
continual change fiom the action of the planets, and therefore 
as deduced fiom observations made m difieient centuries, will 
have different values. 

'J’ho value of p is at present mci easing, or the mean motion 
IB being accelciatcd, and it will continue thus to increase for a 
period of immonse, but not infinite duiation, foi, as shewn by 


♦ Pontccoulant, HxjHibme du Mondo^ yol iv p 600 
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Lagiange, the actions of the planets on the eccentricity of the 
earth’s oihit will be ultimately leversed, e' will cease to 
and begin to mcrease, and consequently p will begin to decrease, 
and tire seculm acceletation will become a seculat retardation 
It IS worthy of romaik that the action of the planets on the 
moon, thus transmitted thiongh the eaith’s orbit, is more con- 
sider able tlian their direct action 


Inequahties dependi7ig on tJte Figure of the Fat th 

100 The earth, not being a perfect sphere, will not attract 
as if the whole of its mass were collected at its centre hence 
some correction must he introduced to take into account this 
want of apheiicity, and some i elation must exist between the 
oblateness and the distmbance it produces Laplace in ex- 
amining Its effect found that it satisfactorily explained the 
mtioduction of a toim in the longitude of the moon, which 
Mayei had discovered by observation, and the argument of 
which IS the tiuc longitude of the moon’s ascending node. 

By a compaiison of the observed and theoretical values of 
the coefficient of this term, we may determine the oblateness 
of the earth with as great, if not greater, accuracy than by 
actual measures on the siuface 

101 By pui suing his investigations, with leference to the 
oblateness, in the expiession for the moon’s latitude, Laplace 
found that it would there give rise to a term in which the 
argument was the true longitude of the moon 

This term, which was imsuspected before, will also serve 
to determine the earth’s oblateness, and the agreement with the 
result of the preceding is almost perfect, giving the compres- 
sion which IS about a mean between the different values 

obtained by othei methods 


* Pontecoulant, SgsUme du Motide^ vol iv 
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Pertm bations due to J’lmi', 

102 After the expicssioji foi the moon’s longitude had boon 
obtained by thcoiy, it was found that tlioic was still a sligJit 
deviation between hci calculated and obscivcd plaeeSj and JBurg, 
who discovered it by a discussion of the obseivatioiis of Lahue, 
Hamstccd, Biadley, and Maskelyue, thought it could be lepre- 
sciitod by an inequality whose peiiod would be 181 yeais aiul 
coefficient 15", This was cntiiely coujectiual, and thougli 
seveial attempts wcie made, it was not accounted foi by thooiy 

About five yeais ago, Piofcssor Hausen, of Secbcig in Gotha, 
having commenced a icvision of the Lunai ''J’licoiy, found two 
teims, which had hithci'to been neglected, due to the action 
of Venus One of them is ducct and aiiscs fiom a ‘leminkable 
^numeiical lelation between the anomabstie motions of the moon 
‘and the sidereal motions of Venus and the caith ; the otlier is 
‘ an indirect effect of an inequality of long pciiod m the motions 
‘of Venus and the caith, which was discovciod some years ago 
‘ by the Astronomer Royal 

The periods of these two inequalities is oxtioinely long, one 
being 273 and the other 230 yeais, and their coefficKuits me 
respectively 27 4" and 23 2". ‘ These are considerable quantities 
‘ in comparison with some of the inequalities already recognised 
‘m the moon’s motion, and, when applied, they are found to 
‘account for the chief, indeed the only lemnining, cmpuual 
‘ poition of the moon’s motion in longitude of any oonsequmico , 
‘so that their discovciy may be consideied as a jiractical c‘om- 
‘pletion of the Lunar Theoiy, at least for tlu' piescuii asUo- 
‘noinieal age, and as establishing the entiic ilominion of the 
‘Newtonian Theory and its analytical application over that 
‘rcfiactory satellite* ’f 


P \ 11 ^^ 848 ^ Annual General Meeting of tlic Ilo> ul Astronomical Society, 

t Address of Sir John Herschel to the Meeting of the Koyal Astronomxcal 
Society 

a 2 



84 


LUNAR THIORY 


Motion of the Ecliptic 

103 We have seen, Art. (14), that our plane of reference 
IS not a fixed plane, but its change ot position is so slow that 
we have been able to neglect it, and it is only when the 
approximation is carried to a high ordei, that the necessity 
arises for taking account of its motion 

It has been found to have an angular velocity, about an 
axis m its own plane, of 48" m a century, and the collection 
thus introduced produces in the latitude of the moon a teim 
— con cos(0 — ^), 

where cd is tlie angular velocity of the ecliptic, - tlie angular 

c 

velocity with which the ascending node of the moon’s orbit 
recedes from the instantaneous axis about which the ecliptic 
rotates, ^ the longitude of this axis at time and 6 the longitude 
of the moon at the same mstant 

Let TAM (fig. 11^) be the position of the ecliptic at time 
A the point about which it is turning, TA = 

NM the moon’s orbit, M the moon, and Mm a perpen- 
dicular to the ecliptic j Tm = 6 , Mm = lat = y9. 
Let APN''m! be the echptic after a time St 
Any point whose longitude is L may be considered as moving 
perpendicularly to the echptic with a velocity g) 8in(L — <j)) 
Hence, if z be the mclination of the oibit, and N the longi- 
tude of the node, the point N will move m the direction NF 
with a velocity to And N* will move along PN' 

with a velocity co cat^, 

dK 

therefore ~di ^ ^ sm(A"~ <f>) cot i. 

Again, the point of the echptic 90“ m advance of will move 
towards the moon’s orbit with a velocity (o sin(90 -f </>) ; 
d% 

therefore - _ o, cos(A'-. <j[>). 

Now, cot f, CD, and - - may be considered constant in 

integrating ; 
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therefore SN = cm coa{N-(f>) cotf, 

S* = cm sin(ii7’— ^), 

where , 

dt c’ 

and if NM = we have 

SN 


Byjt 


cosi 


cm 

sini 


CQa{N—<f)). 


Now, slnyS = sinf smi/r; 

' coSjQ. S)8 = cos* smi/r Sf + sini coa yjr Byjr 

= c6i){cos t sm-i/r sin(iV— <j )) — cos yjr co 3 {N'- 
but cos* sini^ = cos,S am(d~N) and cos’yfr = cos/S C 08 (^— 
therefore Bfi = — cm C08[5 — 

The discovery of this term is due to Piofessor Hansen, its 
coefficient is extremely small, about 1 5" ; but, being of a totally 
different nature from those due to successive approximations, it 
was thought desirable to examine it, and the above investigation, 
which was communicated to me by J. C Adams, Esq , will be 
read with interest on account of its elegance and simplicity 


We may obtam an approximate value of the coefficient cm by 

substituting for it — , where a is the number of seconds thiough 

which the echptic is deflected in one year = 0 48", and n is the 

number of years m which the node of the moon’s orbit makes 

Stt • 

a complete revolution = 18 6 ; for then, — is the angle described 


by the node in one year , therefore, — is the ratio of 


6 ) 


<l>) 
dt ’ 


supposing <j£) to remain constant, which is nearly 


the case , therefore, 


noL 9 3 x 0 48" , ^ 

c© = — = — = 1 42" * 

27r 3 14 


* This affords the solution of a problem proposed in the Senate-House m 
the January Examination of 1852 Question 21, Jan 22 
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NoU on the Nuniencal Values of the Coefficients 


104 When the penods of two of the tenns, m the third 
method given in Art (62), differ but slightly, for instance if 
9 and <f> go through their pei iodic vaiiations very nearly m 
the same time, the method could not then with safety be ap- 
plied , for, smce the same values of 9 and would very nearly 
recur togethei durmg a longer time than that through which 
the observations would extend, the two terms would be so 
blended m the value of V that they would enter nearly as one 
teim — the difference between 9 and ^ would be very nearly 
the same at the end as at the beginning of the senes of obsei- 
vations. 

105 Lot us suppose the penods to be actually identical, 
so that <^ = d + a, a being some constant angle ; then 

B sm 9 + Qsin^ 

may bo wiitten {B+0 cosa) sin0 + (7 sma cos^, 

01 V = A + [B + O cosa) sm^ + C sma cos^ + 

If now we divide the observations, as before, into two sets, 
conespondmg to the positive and negative values of sm0, the 
terms involving cos0 in each set will be as often positive as 
negative, and will disappear in the summation of each set , and, 
foUowmg the process of the metliod, will give 

B+0 cosa = M suppose. 

Dividing again into two sets correspondmg to the positive 
and negative values of cos^, the terms in sind will be cancelled, 
and the same process will give 

0 sma = N suppose. 

Treatmg m the same way with respect to the angle (f>, we 
get two results, 

C + B cosa = M', 

— B sma = N ' , 
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from these four equations we easily get 

„ MN+MN'^^, M'N+MN' 

f7- 

Nj If/ N' are connected by the equation of condition, 
ikP - ir = 


When the periods of 6 and cjE) are nearly, but not exactly, 
the same, this equation of condition will not hold, and the pre- 
ceding values of B and C would not be exactly conect, but 
yet they would be veiy approximate, especially if the mean 
between the two values of B be taken 


106 We may also, after having taken one of these slightly 
eiioneous values for jB, make a further correction by establishing 
as it weie a counterbalancing erroi in the value of G Let 
B' be the value so found for B , then, from the V of each of the 
observations subtiact the value B^ sin0, the icsult U will be 
very nearly equal to A G sin<^ -t- &c , and from the n equations 

— A G + 

ZT^== A G -h 

U,, = A-\-G siiK^,^ + 

a value (7' of G will be obtained, by the rule of Ait (G2), which 
will be veiy approximate, and, at the same time, agiee bcttei 
with B’ in satisfying the equations than G itself would do 

107 When two terms whose periods are nearly equal do 
occm, it IS plain, by examining the values of M and that 
the eriors which would be committed by following the rule, 
without taking account of this peculiarity, would be the taking 
5+ (7 cos a and J5cosa for B and 0 respectively 
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HISTOEY OP THE LUNAE PROBLEM BEFORE NEWTON. 

108 The idea which most probably suggested itself to the 
mmds of those men who first considered the motion of the moon 
among the stars, was that this motion is uniform and circular 
about the earth as a centre 

This first result is represented m our value of the longitude 
by neglectmg all small terms and wiitmg 6 

109. It must, however, have been very soon perceived that 
the actual motion is far firom being so simple, and that the moon 
moves with very different velocities at different times. 

The earliest recorded attempts to take into account the ir- 
regularities of the moon’s motion were made by Hipparchus, 
(140 BC) He imagmed the moon to move with uniform 
velocity m a circle, of which the earth occupied, not the centre, 
but a point nearer to one side. By a similar hypothesis he had 
accounted for the irregularities in the sun’s motion, and his 
success in this led him to apply it also to the moon. 

It 18 dear that, on this supposition, the moon would seem to 
move faster when nearest the earth or m perigee, and slower 
when m apogee, than at any other pomts of her orbit, and thus 
an apparent unequal motion would be produced 

Let BAM (fig, 12) be a circle, OA a radius, B a point in AG 
near 0 j CB, JSB two parallel Imes makmg an angle oc with GA. 

Suppose a body M to describe this circle uniformly with an 
angular velocity the time being reckoned from the instant 
when the body was at B, and the longitude as seen fiom B 
being leckoned from the line BJD, 
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thcrefoie DJEM=6, BOM=pt, 


Now ^.8 
shall have 


AEM=e^a, AOM^p^a. 
a small fraction, and if we represent it by e, we 
EC 

sm -¥= -^miAEM 


= e sin(0 — a) 

“ e sinQ?^H-if- a), 


tan if = 


6Sin(j9^--a) 

1 — e cos(p^ — a) ’ 


this would give if, and then 6 by the formula 6 =pt + M 
This was called an eccentric^ and the value of e was called the 
eocentncity^ which, for the moon, Hipparchus fixed at sin 5® V. 


110 Another method of considering the motion was by 
means of an epicycle^ which led to the same result 

A small circle EM (fig 13), with a radius equal to EC of 
(fig 12), has its centre m the circumfei ence of the circle EFD 
(which has the same radius as that of the eccentric), and moves 
round E with the uniform angular velocity p, the body if being 
carried m the circumference of the smaller circle, the radius PM 
remaining parallel to itself, or, which is the same thmg, revolving 
from the ladius PE with the same angular velocity^, so that 
the angle jEPif equals PEA. 

Now, when the angle AEP equals the angle ACM of the 
former figure, it is easily seen that the two tiiangles jEPif, ECM 
are equal, and theiefore the distance EM and the angle AEM 
will be the same in both, that ,is, the two motions are identical 


111 The value of e being small, we find, rejectmg e\ &c , 
if = e sin(^i — a), 

thcrefoie 0 + e sm(j9i— a). 

If we 1 eject terms of the second order in our expression 
for the longitude, and make c = 1, we get. Art (51), 

0 = 2c sin(^^ — a), 
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which Will he identical with the above if we suppose the eccen- 
tiicity of the eccentric to be double that of the elliptic orbit 

Ptolemy (a d 140) calculated the eccentricity of the moon’s 
01 bit, and found for it the same value as Hippaichus, viz 

sm5° 1' = nearly 

The eccentricity in the elliptic oibit is, we know, about 
These values will pietty neaily reconcile the two values of 6 
given above, and this shews us, that for a few revolutions the 
moon may be considered as moving in an eccentric, and her 
positions in longitude calculated on this supposition will be 
coirect to the first ordei. 

Hei distances from the earth will not however agiee, for 
the ratio of the calculated gieatest and least distances would 

1 _ J i u while that of the time ones would be ^ 

or , which differ by 

It would, thcicfoie, have lequiied two diffcient eccentiics to 
account for the changes in the moon’s longitude and in her 
radius vector Changes in the latter could not, however, be 
easily observed with the rude instruments the ancients pos- 
sessed, and It was very long before this inconsistency was 
detected 

112 We have said that the moon’s longitude, calculated 
on the hypothesis of an eccentiic, will be pietty accuiatc foi 
sl few revolutions 

The data requisite foi this calculation are, the mean angular 
motion of the moon, the position of the apogee, and the mag- 
mtude of the eccentricity' 

But it WHS known to Hippmchus und to the sistionoiners of 
his time, that the point of the moon’s oibit wheie she seems 
to move slowest, is constantly changmg its position among the 
stars Now this point is the apogee of Hippaichus’s eccentiic, 
and he found that he could vciy conveniently' take account of 
this fuither change hy supposing the eccentric itself to have 
an angulai motion about the oaith m the same diiection as 
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the moon herself, so as to make a complete revolution in about 
nine years, or about 3° m each revolution * 

This motion of the apsidal line follows also from our expres- 
sion for the longitude, as shewn m Art (66) It is theie, 
however, connected with an ellipse instead of an eccentric ; and 
though the discovery that the elliptic is the true foim of the 
fundamental orbit was not the next in the ordei of time aftei 
those of Hipparchus, yet, as all the irregularities which were 
discoveied m the intervening seventeen centunes are common 
both to Hipparchus’s eccentric and to Kepler’s ellipse, it will 
be as well for us to consider at once this new form of the 
orbit 


JElhptio Fo'tm of the 0?htt 

113. We need not dwell on the steps which led to this 
great and important discovery Keplei, finding that the pie- 
dicted places of the planet Mars, as given by the circular 
theories then m use, did not always agiee with the computed 
ones, sought to reconcile these variances by other combmations 
of circular orbits, and after a great number of attempts and 
failures, and eight yeais of patient mvestigation, he found it 
necessary to discard the eccentrics and epicycles altogethei, 
and to adopt some new supposition An elhpse with the sun 
m the focus was at last his fortunate hypothesis, which was 
found to give lesults in accoi dance with observation, and this 
form of the oibit was, with equal success, afterwaids extended 
to the moon as a substitute foi the eccentric . but the departures 
from elliptic motion, due to the disturbing force of the sun, 
are, m the case of the moon, much greater than the distuib- 
ances of the planet Mars by the other planets 


* On the supposition of an epicycle, this motion of the apse could as easily 
be lepresented by supposing the ladms which connects the moon with the 
centre of the epicycle to have this unifoim. angulai velocity of about 3® m each 
1 evolution, and also in the same direction 
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In Kepler’s hypothesis, then, the earth is to be considered 
as occupying the focus of an ellipse, in the perimeter of which 
the moon is movmg, no longer with either uniform hnear or 
angular velocity, but in such a manner that the radius vector 
sweeps over equal areas in equal times 

This agrees with our investigation of the motion of two 
bodies. Art (10) 


JBvection 

114. The hypothesis of an eccentric, whose apse line has 
a progressive motion, as conceived by Hipparchus, served to 
calculate with considerable accuracy the circumstances of 
eclipses , and observations of eclipses, requirmg no instruments, 
were then the only ones which could be made with sufficient 
exactness to test the truth or fallacy of the supposition 

Ptolemy (A D 140) having constructed an mstrument, by 
means of which the positions of tlie moon could be observed m 
other parts of her orbit, found that they sometimes agreed, but 
were more frequently at variance with the calculated places, 
the greatest amount of error always taking place at quadrature 
and vanishing altogether at syzygy. 

What must, however, have been a source of great perplexity 
to Ptolemy, when he attempted to investigate the law of this 
new irregularity, was to find that it did not return m every 
quadrature, — in some quadratures it totally disappeared, and 
in others amounted to 2 ° 39', which was its maximum value 
By dmt of careful comparison of observations, he found that 
the value of this second inequality m quadrature was always pro- 
portional to that of the first in the same place, and was additive 
or subtractive accoiding as the first was so and thus, when 
the first inequality in quadrature was at its maximum or 5"* 1', 
the second increased it to 7 ° 40', which was the case when the 
apse line happened to be in syzygy at the same time* 

♦ It would seem as if Hipparchus had felt the necessity for some further 
modification of his first hypothesis, though he was unable to determme it, for 
tHore IS an observation made by him on the moon m the position here specified 
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But if the apse line was in quadrature at the same time as 
the moouj the second inequality vanished as well as the first. 

The mean value of the two mequalities combined was there- 
fore fixed at 6“ 20^'. 

115 To repiesent this new mequality, which was subse- 
quently called the Evection^ Ptolemy imagined an eccentric in 
the circmnference of which the centre of an epicycle moved 
while the moon moved in the circumference of the epicycle. 

The centre of the eccentric and of the epicycle he supposed 
in syzygy at the same time, and both on the same side of 
the earth 

Thus, if E repiesent the earth (fig 14), 

8 . . sun, 

M moon, 

c the centre of the eccentric RKT in syzygy, 
jB, the centre of the epicycle, would also be in syzygy. 

Now conceive c, the centre of the eccentric, to describe a 
small circle about m a retrograde direction cc', while i?, the 
centre of the epicycle, moves m the opposite direction, m such 
a manner that each of the angles 8Ec^ 8ER' may be equal 
to the synodical motion of the moon, that is, her mean angular 
motion from the sun; 8E8* bemg the motion of the sun in 
the same tune. 

Now we have seen. Art (110), that the first mequality was 
accounted for by supposing the epicycle RM to move mto the 
position rm, r and R being at the same distance from jS, and 
rm parallel to RM^"^ the first inequality being the angle rEm 
But when the centre of the epicycle is at iJ', and EM’ is parallel 
to the mequality becomes R’EM'^ and we have a second 
correction or mequality mEM’. 

"wlien the error of his tables 'would be greatest , and at a tune also -when she 
was m the nonagesimal, so that any error of longitude, arising from her yet 
uncertain parallax, would be avoided Ptolemy, who records the observation, 
employs it to calculate the evection, and obtains a result agreeing with that 
of his own observations (See Dmambre, Ast A7icim7ie ) 

* For simphcity we leave out of consideration the motion of the apse 
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116 That this hjpothesis will account for the phenomena 
observed hy Ptolemy, Art (114), will he readily undeistood 
At syzygies, whether conjunction or opposition, the centies 
of the eccentric and epicycle aie m one hue with the eaith 
and on the same side of it , the points i and E comcide, as also 
m and M. Hence mEM = 0 

At quadiatures (figs 15 and 16) d and E are m a straight 
line on opposite sides of the earth, and therefoie E and r at 
their furthest distance If, however, M.' and m be at the same 
tune m this line, or, m other woids, if the apse line be in 
quadratmes (fig 15), the angle mEM' wiU still be zero, or there 
will be no erioi m the longitude But, if the apse line is in 
syzygy (fig 16), the angle mEM' attams its greatest value* 
Ptolemy, as we have said, found this greatest value to be 
2° 39', the angle mEt being then 5° 1' 


117. Copernicus (i. d 1543), havmg seen that Ptolemy’s 
hypothesis gave distances totally at variance with the obser- 
vations on the changes of apparent diameter ,t made another 
and a simpler one which accounted equally well for the m- 
equahty m longitude, and was at the same time more correct 
m its representation of the dia ta nc e s 

Let E (fig 17) be the earth, OD an epicycle whose centre 
G describes the circle COO" about E with the moon’s mean 
angular velocity 

Let 00, & ladius of this epicycle, be parallel to the apse 
line EA, and about 0 as centre let a second small epicycle be 
described, the radii GO and OJWhemg so taken that 


GO-OM , CO+OM 

— G^=sm5 1,and— 


= sm7®40\ 


♦ H Ptolemy tad used the hypothesis of an eecentnc instead of an epicycle 
for the first inequality of the moon, an epicycle Tvould have represented the 
secoad inequality more simply than his method did Dr Whewell’s ffutory 
of the Inductive Sciences, toI i p 230 

t See Delamhre, Ast Modeme, yol i* p 116 "Whewell’s History of 
Inductive Sciences, toI i p 395 
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The ladius OM must now be made to revolve from the radius 
OG twice as lapidly as ^(7 moves from JES^ so that the angle 
COMmeiy be always double of the angle CU8. 

From this construction, it follows that m syzygies the angle 
CJS8 being 0° or 180°, the angle COM is 0° or 360°; and there- 
fore G and M are at their nearest distances, as in the positions 
C' and G'” m the figure Then G3€ = CO — OM^ and the 
angle GEM will range between 0° and 5° 1', the greatest value 
bemg attamed when the apse Ime is m quadiatuie 

When the moon is in quadiatuie GE8 = 90° or 270°, and 
theiefore, COM— 180° oi 540° and G and M aie at their gieatest 
distance apart, as in the position (7" , then, CM = (70 -1- OM^ 
and the angle GEM will range between 0° and 7° 40', the foimer 
value when the apse line is itself m quadrature, and the latter 
when it IS in syzygy 


118 Thus the lesults attained by Ptolemy’s construction 
are, as far as the longitudes at syzygies and quadiatures aie 
concerned, as well lepresented by that of Copeimcus, and the 
variations m the distances of the moon will be far more exact, 
the least apparent diameter being 28' 45" and the gieatest 37' 33"; 
whereas, Ptolemy’s would make the gieatest diameter 1° ^ 

The values which modem obseivations give vaiy between 
28' 48" and 33' 32" 

119 It will not now be difiicult to shew that the introduc- 
tion of this small epicycle conesponds with that of the term 
^me sm{(2 - 2??i — c);pt- 2/3 + a} in oiu value of 6 

Foi, lefeirmg to (fig 17), we have 


OEM = sin OEM = ^ sm OME 
= ^sa^{GOM-AEM) 

= 8EG-AEM) 

— sin{2(moon’s mean long — sun’s long ) 
- (moon’s''^r^^e long. — long of apse)}. 


* Delambie, Ast Modeme 
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and OEM being a small angle whose maximum is 1® 19^', we 
may write moon’s mean longitude mstead of the true m the 
argument, and also EG for OE^ therefore, 

OEM — sm[2(moon’s mean longitude — sun’s longitude) 

— (moon’s mean longitude — longitude of apse)] 


= ^)] — {pt— (1 — c)^^ + a}] 

= 4770" sm{(2 — 2m— c)jpt— 2;3 + a}. 

The value of the coefficient is from modern observations found 
to be 4589 61" 


120 In Art (70), we have considered the effect of this 
second inequality in another light, not simply as a small quantity 
additional to the first or elliptic mequality, but as formmg a part 
of this first; and theiefore, modifying and constantly altering 
the eccentricity and the uniform progression of the apse line. 

Boulliaud (a d. 1645), by whom the term Evection was first 
applied to the second inequality, seems to hint at something 
of this kind m the rather obscure explanations of Ms lunar 
hypothesis, which never having been accepted, it would be 
useless to give an account of.* 

In Ptolemy’s theory. Art (115), the evection was the result 
of an apjparent mcrease of the first lunar epicycle caused by 
its approaching the earth at quadratures; but, in this second 
method, it is the result of an actual change in the elements of 
the elliptic orbit. 

D’Arzachel, an Arabian astronomer, who observed m Spam 
about the year 1080, seems to have discovered the unequal 
motion of the apsides, but his discovery must have been lost 


♦ Aprfes avoir ^tabL. les mouvemens et les fipoques de la lime, Boulliaud 
revient k Texplication de r^vection ou de la seconde in§galit6. Si sa tbSorie 
n'a pas fait fortune, le nom du moms est rest6 ‘ En mime temps q.ue la luue 
‘ avance sur son. c6ne autoux de la terre, tout le sysUme de la lune est d&placi ; 
‘la terre emportant la lune, re 3 ette lorn d*eUe Tapogle, et rapprocbe d'autant 
‘ le perigee , mais cette Ivection k <ieB homes fixles ' 

Delambre, JSist de VAst Mod , tom u p 157 
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Sight of, foi Horrocksj about 1640, rediscovered it ‘m conse- 
^ quence of his attentive observations of the lunar diameter he 
^ found that when the distance of the sun from the moon’s apogee 
‘ was about 45“ or 225°, the apogee was more advanced by 25° 

^ than when that distance was about 135° or 315° The apsides, 

‘ therefore, of the moon’s orbit were sometimes progressive and 
‘ sometimes regressive, and reqmred an equation of 12° 30', some- 
^ times additive to their mean place and sometimes subtractive 
‘ from it 

Horrocks also made the eccentricity variable between the 
limits 06686 and 04362. 

The combmation of these two suppositions was a means of 
avoiding the introduction of Ptolemy’s eccentric or the second 
epicycle of Copermcus . their joint effect constitutes the evection 


Variation 

121. After the discovery of the evection by Ptolemy, a 
period of fourteen centuries elapsed before any further addition 
was made to our knowledge of the moon’s motions Hipparchus’s 
hypothesis was found sufficient for eclipses, and when corrected 
by Ptolemy’s discovery, the agreement between the calculated 
and observed places was found to extend also to quadratures; 
any slight discrepancy bemg attributed to errors of observation 
or to the imperfection of instruments 

But when Tycho Brah6 (a d 1580) with superior mstiiments 
extended the range of his observations to all intermediate pomts, 
he found that another mequahty manifested itself Having 
computed the places of the moon for different parts of her orbit 
and compared them with observation, he perceived that she was 
always m advance of her computed place from syzygy to quad- 
lature, and behmd it from quadrature to syzygy, the maximum 
of this variation taking place m the octants, that is, m the 
points equally distant from syzygy and quadrature The moon’s 
velocity therefore, so far as this mequahty was concerned, was 


* Small’s Astronomical Discover tes of Kepler , p 307 
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greatest at new and full moon, and least at the fiist and thud 
quarter * 

Tjcho fixed the maximum of this inequality at 40' 30" The 
value which lesults fiom modem observations is 39' 30" 

122 We have already two epicycles, or one epicycle and an 
eccentiic, to explain the fiist two inequalities* by the inti educ- 
tion of anothei epicycle or eccentric, the variation also might 
have been bi ought into the system, but Tycho adopted a dif- 
ferent method f hhe Ptolemy, he employed an eccentiic foi the 
cvection, but for the first or elliptic inequality he employed a 
couple of epicycles, and this complicated combmation, which it 
is needless further to desciibe, lepresented the change of distance 
bettei than Ptolemy’s 

To introduce the mrvxtwn,^ he imagmed the ccntie of the 
laigei epicycle to librate backwaids and forwaids on the ecccn- 
tiic, to an extent of 40|' on each side of its mean position, this 
mean place itself advancmg unifoi’mly along the eccentiic with 
the moon’s mean motion in anomaly , and the libi ation was so 
adjusted, tliat the moon was in her mean place at syzygy and 
quadrature, and at hei furthest distance fiom it in the octants, 
the period of a complete libration bemg half a synodical i evo- 
lution 


* < It appears that Mohammed-Aboul Weih al-Bouzdjani, an Arabian astro- 
‘ nomer of the tenth century, who resided at Cairo, and observed at Bagdad m 
‘ 976, discovered a third inequality of the mpon, in addition to the two cx- 

* pounded by Ptolemy, the equation of the centre and the evection This 
‘ third mequahty, the variation, is usually supposed to have been discovered 
‘ by Tycho Brahe, six centuries later In an almagest of Aboul-Wcfa, a 
‘part of which exists in the Royal Library at Pans, after describing the two 
‘ mequahties of the moon, he has a Section IX , “ Of the third anomaly of the 
‘moon called MvJumaX or Prosnmsu" But this discovery of Aboul-Wefa 
‘ appears to have excited no notice among his contemporaries and followers , at 
‘ least it had been long quite forgotten, when Tycho Brah6 rediscovered the 

* same lunar inequality Whewell’s Hist of Isiductive ScienoeSf vol i p 243, 

t For a full description of Tycho’s hypothesis, see Delambre, Bist de VAst. 
Jlfod,tom I p 162, and An Aeeo«>t<q^<AcA«i)'o»om8ca7X)i«cormeso/A«t)fcr, bv 
Robert Small, p 135 
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Annual EqtiaUon 

123 Tycho Biahd was also the disco veiei of the fomth 
inequality, called the annual equation This was connected 
with the anomalistic motion of the sun, and did not, like the 
piovious mequalities, depend on the position of the moon m hei 
orhit 

Havmg calculated the position of the moon conesponding to 
any given time, he found that the ohseived place was behind 
her computed one while the sun moved from peiigee to apogee, 
and befoie it m the othei half year 

Tycho did not state this distinctly, but he made a conection 
which, though wiong m quantity and applied m an indiiect 
mannei, showed that he had seen the necessity and understood 
the law of this mequality 

lie did not try to represent it by any new eccentric oi 
epicycle, but ho mci eased by (8m 13s,) sin(.swM’9 anomalif) the 
tunc which had seivod to calculate the moon’s place thus 
assuming that the true place, after that interval, would agice 
with the calculated one Now, as the moon moves tliiough 4' 30" 
111 8m 13s , it IS clcai that adding (8m 13s ) sm(si^ 72 ’s anomaly) 
to the tune is the same thing as subti acting (4' 30") sm(sM9i’5 
anomaly) fiom the calculated longitude, which was theiefoie 
the collection virtually intioduced by Tycho f Modern obser- 
vations shew the coefficient to be 11' 9" 

Wo have seen, Ait (75), how this inequality may be infciicd 
fiom our equations 


Eeductwn 

124 Tho next mcquahty in longitude whidi we have to 
consider, is not an inequality in the sanio sense as tho foregoing , 

* That IS, tho eq.uatiou of time which he used loi the moon diffoied by that 
quantity from that used foi tho sun 

t Horrocks ( 1639) made the correction in the some manner as Tj cho,jbut so 
increased it that the corrcsiiondmg coefficient was 11' 61" instead 
Flamsteed wos the flist to apply the collection to the longitude instead of tho 
time 



100 


LUNAR THEORY 


that 1S5 it does not aiise from any irregularity in the motion of the 
moon herself m hei orbit, but simply because that orbit is not 
in the same plane as that in which the longitudes aie icckoncd, 
so that even a regular motion m the one would be neccssaiily 
inegulai when lefened to the othei Thus if NMn (fig 10) 
be the moon’s oibit and TNm the ecliptic, and if M the moon 
be leferred to the ecliptic by the great circle Mm pci pondicular 
to it, then ICVand mN are 0®, 90°, 180°, 270°, and 360° simul- 
taneously, but they diiSfer for all intermediate values the 
difference between them is called the reduction 

The difference between the longitude of the node and that 
of the moon m her orbit being known, that is the side NM of 
the light-angled spherical tiiangle iVTIfm, and also the angle 
N the mclination of the two mbits, the side Nm may be calcu- 
lated by the rules of spherical tiigonomctiy, and the difference 
between it and applied with a piopei sign to the longitude 
in the 01 bit, gives the longitude m the ecliptic 

Tycho was the fiist to make a table of the 1 eduction instead 
of calculating the spherical triangle His foimula was 

reduction = tan"|-/ sin2X - -^tan^^/sin iX, 

where I is the inclination of the 01 bit and L tlie longitude of 
the moon dimmished by that of the node 

The first term corresponds with the term — ^/r mu2{(;pf — y) 
of the expiession for 6 


Latitude of the Moon 

125 Ttat the moon’s orbit is inclined to tlic ecliptic was 
known to the eailiest astionomcrs, fiora the the non-rcciinoncc 
of eclipses at every new and full moon; and it was also known, 
smee the eclipses did not always take place m the same pait 
of the heavens, that the Ime of nodes icprescntcd by Nn 
(% 10) has a letrograde motion on the ecliptic, N moving 
towaids T 

Hippaichiis fixed the inclination of the moon’s mbit to tlic 
ecliptic at 5°, which value he obtained by obscivmg the gicatost 
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distance at which she passes to the north or south of some star 
known to be m or very near the ecliptic, as for mstance the 
bright star Eegulus; and by comparing the recorded eclipses 
from the times of the Chaldean astronomers down to his own, 
he found that the line of nodes goes round the ecliptic m a 
retiograde direction m about 18| years 

This lesult is indicated in our expression for the value of the 
latitude by the term h sin(^0 “7)? ^.s we have shewn Art (78) 

126 Tycho Brah^ further discovered that the inclination of 
the lunar oibit to the ecliptic was not a constant quantity of 5“ 
as Hipparchus had supposed, but that it had a mean value of 
5° 8', and ranged through 9' 30" on each side of this, the least 
inclination 4° 58^ occurimg when the node was in syzygy, and 
the greatest 5° 17^ being attained when the node was m 
quadiature * 

He also found that the retrograde motion of the node was 
not uniform* the mean and true position of the node agreed 
very well when they were m syzygy or quadratuie, but they 
weie 1° 46' apart in the octants 

By referring to Art (80), we shall see that these collec- 
tions, introduced by Tycho Biah^, correspond to the second 
tcim of our expiession for s 

Since Hipparchus could observe the moon with accuracy 
only in the eclipses, at which time the node is m or neai syzygy, 


*** Ebn JoTinis, an Ai'i'bian astronomer (died a d 1008), whose works 
were translated about 30 >ears since by Mens Sedillot, states that the in- 
clination of the moon’s orbit had been often observed by Aboul Hassan Aly- 
ben-Aonajour about the year 918, and that the results he had obtamed were 
generally greater than the 5® of Hipparchus, but that they mned comxderaUy 
Ebn Joums adds, however, that ho himself had observed the inclination 
several times and found it 5® 3', which leads us to infer tliat he always observed 
in similar circumstances, for otherwise a variation of nearly 23' could scarcely 
have escaped him See Delambre, Hist de VAbt du Moyen Age^ p* 139 

The mean value of the inclination is 5® 8' 55 46", —the extreme values are 
4® o7' 22' and 5" 20' 6" 

The mean daily motion of the line of nodes is 3' 10 61", oi one revolution 
in 0793 30 days, 01 18y 218d 21h 22m 46 b 
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we see why the value he found for the inclination of the orbit 
was appioximately its minunum value, and also why he waa 
unable to detect the want of uniformity m the motion of tho 
node, 

127 To represent these changes in the position of tho 
moon’s orbit, Tycho made the following hypothesis. 

Let ENF (fig 18) be the ecliptic, K its pole, BAO a small 
circle, havmg also K for pole and at a distance fiom it equal 
to 5° 8'. Then, if we suppose the pole of the moon’s oibit to 
move uniformly in the small circle and in the direction BAG, 
the node N, winch is at 90° from both A and K, will retiogiade 
imifoimly on the ecliptic, and the inclination of the two orbits 
will be constant and equal to AK 

But mstead of supposmg the pole of the moon’s orbit to be 
at A, let a small ciicle aSccZ be described with A as pole and 
a radius of 9' 30", and suppose the pole of the moon’s oibit to 
describe this small circle with double the velocity of the node 
m Its synodical revolution which is accomplished in about 
346 days, in such a manner that when the node is in syzygy 
the pole may be at a, the nearest point to JC, and at c the 
most distant point when the node comes to quadrature, at 5 in 
the first and third octants, and at d in the second and fourth, 
so as to descnbe the small cncle in about 173 days, the centre 
A of the small circle retrograding meanwhile with its uniform 
motion 

.?y this method of representmg the motion, we see that 
when node is m syzygylthe indmationjZa = 5° 8' - 9^' = 4° 58^', 
quadratmej of the oibit is [Kc = 5° 8' + 9J^' = 5° 17 ^'^ 
while at the octants it has its mean value Kb = Kd = KA = 5° 8' 

Agam, with respect to the motion of the node, since N is 
the pole of KaAc, it follows that when m syzygy and quad- 
rature, the node occupied its mean place, in the fiist and third 
octants, the pole bemg at b, the node was behind its mean 
place by the angle bKA = (9' 30") cosec .5° 8' = r46', neaily, 
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and it was as much m advance of its true place in the second 
and fomth octants 

So that the whole motion of the node, and the correction 
which Tycho had discovered, were properly represented by 
this hypothesis, which is exactly similar to that which Copei- 
nicus had imagined to explain the precession of the Equinox. 


THE END 
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fore, heieafter omit such terms as -4 cos (3 — B) altogether, and 
merely write 

M = a(l + e cos(c0 — a)} + 

38 So fai, all that we know about c is that it differs from 
unity at most by a quantity of the hrst order, but its value wiU 
be more and more conectly obtamed by always writing, m the 
successive appioximations, a + ae cos {c6 — a) for the first two 
terms of the value of w, then the coefficient of cos(c0 — a) m 
the differential equation must equal ae(l — c^); and this will 
enable us to determine c to the same order of approximation 
as that of the differential equation itself See Arts (48) and (94) 

39 In cai lying on the solution of s, the same difficulty 
anses as in w, and it will be found necessary to change it into 

5 « i &m{g9 - 7) . , 8\^ 

ff being a quantity which diffeis from unity at most by a quan- 
tity of the fiist order See Aits (49) and (95) 

40 The equation will also be modified by this change 
in the value of w, 

dt ^ I 1 

dd ” ha^ {1 4 e cos(c0— a)}'* ’ 

det 1 1 

deff"^ ha^ {1 4 ecos(c6-“ a)}** 

Here c9 and ct hold the places which 9 and t occupied 
in (33) , therefoie 

c9 cpt + 2e Bm(cpi-- a), 

01 ^ 4 sin(qp^- a) . 

to the first Older, since - = e to the first order 

41. Smee the disturbing forces are to be taken into account 
in the next approxunation, we shall have to use the value of 
u' found m (18), which is 

w' = a {1 4 d cos(fl' — §)} 
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but this mtioduces d\ we must therefore furtliei modify it by 
substituting for 6' its value m terms of and it will be found 

sufScient, for the purpose of the present woik, to obtain the 

connexion between them to the first orderj which may be done 
as follows 

Let 7n be the ratio of the mean motions of the sun and moon, 
p their mean angular velocities , p' = wj?, 

4 pt mean longitudes at time l3 being the sun’s 

longitude when if = 0, 

6 true longitudes at time tj 
^5 a longitude of perigees when ^ = 0 , 

thciefoie 6' — f = sun’s tiue anomaly, 

and f z= mean anomaly 

But, by Art (13), 

time anomaly = mean anomaly + 2e' sin (mean anomaly) + &c , 
theiefore ff = p't + ^ + 2e' sin 4 )8 - ^ 4 

= 7npt 4/8 + 2e sin {m 2 )t 4 /S - 4- 

= md 4 ^ 4 2e' sin {tnd 4 /8 — f ) 

to the fiist Older, 

because 0 — 2e sin {c0 — a) to the first order by (40) 

Whence = a' {1 4 e' cos (w6^ 4/8— f)} to the fiist Older 

42 The values of sm2(0- 6') and cos 2(0-0') can also be 
readily obtained to the same older 

sin 2(0 — 0') 

= sin {(2 - 2m) 0 — 2/3 — 4e' sm(m0 4 /8 — f )} 

= sin {(2 - 2m) 0 - 2/8} — 4e' sin(m04 /S - f) cos {(2 — 2m) 0 — 2/3} 
= sin {(2 - 2m) 0 — 2/8} — 2c' sin {(2 — m) 0 — /8 — f} 

4 2e' sin {(2 - Bm) 0-3/84 5'} 
Similarly, cos 2(0 - 0') — cos {(2 — 2m) 6 — 2/3} 

-2c' cos{(2-7w) 0-/3— f} 4 2e' cos {(2- 3m) 0-3yS-i-5’} 
The first teim of each of these is all we shall require 

d2 
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Section III 

To solve the Equations to the Second Order 

43 Let us recapitulate the results of the last appioxmiation 
[1 + 6 cos(c5 - a)}, 
v! = a {1 + e Q.o^[md + ^ — S’)}, 

5 = i 8111(^0 — 7)3 

0 — 0' = (1 — m) 0 — ^ — 2e' sm(9n0 + )8 — f) 

These values must now be substituted m the expiessions foi 
P T Ps-S T du T ds (d^u \ {T . 
AV’ AV’ AV ’ AVrf0’ AVrf0’ \dd^^'^) ^ 

( d^8 "N r 2^ 

^ + A \j^ ^^3 letammg terms above the second order, when, 

according to the ciiteiion of Art (29), they promise to become 
of tlie second older aftei integrating 

The equations (/S') and (7 ) of Art (30) will then assume 
the forms 

TTZ/IV 

and the integration of these will enable us to obtain u and s 
to the second order , aftei which, equation (a) of Art. (20) will 
give the connexion between 0 and t to the same oidei, 

44 The quantity , which we shall meet with as a coef- 
ficient of the terms duo to the disturbing foice, can be replaced 
by m“a, m bemg the ratio of the moan motions of the sun 
and moon; for, as in Ait (24), 

wi'a'** __ _ (periodic time of moon about earth)'* 

AV fJba^ (periodic time of sun about earth) ^ 

nearly. 
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and 


P' _ time of moon about eaith , 

” pel iodic time of sun about earth ’ ^ ’ 


tlieiefoie 


may be leplaced by m^a ^ 


45 We have thercfoie 


P 


- ) 


m u 

w 


{j +^cos2(0-e')} 


= a [I — sin‘'^(^<70 — 7)} 


cos(m^+y3"- S')in 
AV (l+e cos(c0— a)}^ J 

[| + |oos{(2 — 2m) 6 — 2/3]] 


r 1 - + |/c‘" cos 2(,9'0 — 7) 

= a\ - {1 + 3e cos(mi9-hy8~- ?)) {1 - 3e cos(c0 - a)} 

I [1 + 3 cos [(2 — 2m) 0 — 2^)] 

fl~ f + iJc^ cos 2 (,g^0~7) - [1+3 cos{(2-2?}i) 0 - 2^}] 

= a(— pnV cos [m9 + /3 — f) + -^-mV cos — a) 

[ + \rn^e cos((2 — 2/>^ — c) 0 — 2/3 + a] 


^ It may perhaps be imagined, smce the 01 bits described arc not acciuatcly 
ellipses, and, even it they were, a' and a would be the reciprocals of the scmi- 

, _ - mV* , 

lateia recta, and not those of the bcnii major a\cs, tliat thcicfoio is only 

approximately equal to , and that if it were required to carry our inves- 
tigations to a higher oidei than we propose to do, it would be necessary to 
correct this equation 

But the very fact that the moon’s orbit is not an exact ellipse, shews that 
the yet unknown quantity a is not a definite magnitude whose value is fixed 
at the outset like those of m' and /u It is a constant, it is true , but it is 
one whose value may bo whatever we please, provided the assumption do not 
iiiteifeio with our continuous appioximation 

Now from Alt ( 30 ), where « was first introduced, we sec that it is equal 


to ~ . A* being itself an arhitimy constant brought m by the integration, and 

this bemg the only relation between the quantities 7i and a, wo arc at liberty 

to establish any second one The relation is this second assumed 

relation The reasoning in the text docs not piove this i elation, since, as -v^e 
say, it IS an aibitiary one, but it suggests it The actual value of a is now 
fixed, and the method ol obtaining it from observation ^\lll be cxi)laiiied m tlie 
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The last three terms are retained, though of the thud oidei, 
accordmg to Art (29) * The first of the throe will not rise in 
importance m the value of w, hut it is retained for its subsequent 
use m findmg when it will become of the second oidei The 
other terms of the third ordei, which aiise m the development 
of the expression, are neglected, as the coefficients of Q iii then 
arguments are neither small nor near unity 


T ^ 

AV 2A^w^ 


sm2(0— &) 


__ Z'fYicT {l + cos (m6 4“ /3 — ^)Y 


sin((2 — 2m) 6 - 2/3} 


— - J 


27iV {1 + e cos {c6 — a)]'^ 

== — |m^ {1 + Be' cos(m0 + yS — f)} 

{1 — 4e co3(c0 “■ a) + lOe'* coaled — a)} sin{(2 — 27?z) 6 — 2yS} 
= — f {1 + 3e' cos(m0 + /3 — f)} 

[sin{(2 - 2m) 6 — 2^5} — 2e sin{(2 — 2m — c) 0 — 2/3 + a} 

+ fe'* sm[(2 — 27?i ~ 2c) ^ - 2/3 + 2a}] 
'sm[(2— 2m) 6 - 2/3} - 2c sm{(2 - 2m — c) 0 — 2/3+ a} 
,+ sin{(2 — 27?z — 2c) 0 — 2)8 + 2a} 

We have, in the couise of the reduction, dropped those terms 
which, according to Ait (29), could not produce important terms 
in the resulting value either of w or of ^ The last teim, though 
of the fourth order, is retained because 2 — 2m — 2c is small 

^ ~ (previous expression) {— asc sin(c0— a)} 

= %rr^ae cos[(2 - 2m — c) 5 — 2)8 + a}, 
foi to the first order c — 1 ; 


next chapter If we had assumed = m*Xo, where X is any given quantity 

differing very shghtly from umty, we could still have proceeded with oui 
approximations , for the corresponding value of a, though slightly different 
from the former, would stiU, when used in the equation Z7^', give a value 
of u approximate to the first order 

But if X differed considerably from unity, the value of a which this equa- 
tion would furnish, would no longer render an approximation to the first 

order, and if we made use of the equation we should have no guarantee 

that the next step would he a closer appioximation 
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T ds 

55 ~ (same expiession) [l^g cos [gB- 7)} 

= — sm{(2 — 2 ot — 9) 5 — 2^ + 7}, 

for 7 = 1 to the first ordei 


but 


2 -2m ^ 

Op 

( — - — cos{(2 - 2?» — c) 0 — 2/3 + a} 

+ 2 ( 2 - 2 ! -2c) ^ -2m- 2c) B-2fi + 2a] . 

^ = + + teims of first order, 

2 -2m ^ ’ 


— r — r = — 2e + terms of second order, 

2 — 2m-c ’ 

5e^ ^ 5e^ 

2 (2 — 2m — 2c) 4 (1 — m — c) 

Here the denominator is of the first ordei, and cannot be 
further simplified without a more accurate knowledge of the 
value of c We shall find in the next value of w, Ait (48), 
that 1 — c IS of the second order, and as this result is obtamed 
independently of the teim wo are here Considering, which is 
only retained for the sake of finding there is no impropriety 
in anticipatmg thus far m order to simplify this coefficient, 
which then becomes 

— ^ + terms of the second order , 

4m ' 

• J = fm“ cos{(2 — 2m) 9 - 2^] 

— 3m“<3 cos{(2 — 2m — c) 9 -^213 + a] 

— ^me^ cos{(2 — 2m — 2o) 0 — 2)8 + 2a}, 

Also, by Art (36), 

^ 4- w = a -f quantities of the second order, 
d^s 

+ s = small quantity of the second oidoi at least , 
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— Gm^ae cos {(2 — 2vi ~ c) 0 - 2^ -{- a} 

— ^me~a cos {(2 — 2m — 2c)6 — 2^ -}- 2a}, 

^ (^ "*” ^) / ~ ordei 


Lastly, 

Ps-- 8 


m’su'^ 

“AV" 


(t + f cos 2(0-0')} 


= — [1 -{“ cos {(2 — 2m) 0 — 2/3|] 

{1 + Zg cos (m0 + /3 — 5') “ cos (c0 — a)} 

= - ^m'^h [sm( 9^0 — 7) - "I sm {(2 — 2m — 9 ) 0 — 2/3 + 7 }] 

In all these expressions we have rejected those teims of the 
thud and highei oiders which, according to Art (29), would 
not influence tlie second older ^ 


46 We must now substitute these values m the diffeiential 
equations foi u and 5 , and then integrate, omittmg the com- 
plementary term A cos (0 — 5) , for though, by the theory ot 
ditfeiential equations, this would fonn a necessaiy part of the 
solution, it cannot, as we have seen Art (35), be a pait of the 
coirect values of 01 s 


* Instead of the forces which really act on the moon, wc originally sub- 
stituted three equivalent ones, J% /*, S, those agtun aie, by the picccding 
expressions, replaced by a set of others For, wc miy conceive each of the 

p 

teims m &c to correspond to a foicc, — a component of P, T, oi S , each 

force having the same argument as the teim to which it coriesponds, and there- 
fore going through its cycle of values in the same time Now, by Art (29), 
when the coefficient of 6 in the argument is neai unity, the term becomes 
important m the radius vector, and when near zero, in the longitude hence, 
a force whose period is nearly the same as that of the moon, produces im- 
portant effects, m the radius vector, and a foice whose period is vciy long 
will be important in its effect on the longitude 

See Aiuj\ Tiacts, Planctduj Theory, p 7S 
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47 Since the foim of the solution is known, the actual ex- 
pressions tor u and $ will be obtained with more facility by 
assuming them with arbitiary coefficients, the values of which 
aie aftciwaids detei mined by substitution 

Wo must remember, however, that the coefficients of 
cos(c0 — a) m u and of sin (^0—7) m s must be assumed the 
same as in the fiist approximate solutions, and that these 
assumptions will enable us to obtain the values of c and g to 
the same oidcr of approximation as that to which we are 
woiking, Alt (38) 

48 Considering, firstly, the equation m w, we have 

d^u _ P ^ ^ _ 9 ^ r ^ ,7ft 

<70^ “ “ W “ W ^ “ j J W 

/I 

-h Im^e cos(c0 — a) 

-b lU cos 2 [gO - 7) 

= a ( “ cos {(2 — 2w) 0 — 2/3} 

-b cos{(2 — 27?i - c) 0 - 2)S + a] 

- ■|77iV cos(m0 + iS— ?) 

^-b '^me^ cos {(2 — 2m - 2o) 0 — 2)8 -h 2a} 

The last two terms would not be retained if we wished to 
fin d the value of u only, but, in finding t afterwards, they will 
become of the second older 

^ 1 — f /s'* “ 

+ e cos*(c0 — a) 

+ A cos2(^0““7) 

Assume 71 = a{ + B cos {(2 - 2m) 6 — 2/3} 

-1- C cos{(2 — 2m - c) 0 — 2/3 -b a} 

+ D cos(m0-b/3- ?) 

-b E cos {(2 — 2771 — 20) 0 - 2)8 -b 2a}, 
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Then, by substitution, 

e (1 — c'^) = 

jS {1 — (2 — 2mY] — — 

<7 {1 -- (2 — 2m — c)"*} = 

Z> (1 — m^) = - 

jE {1 — (2 — 2m — 20)*"*} = ; 

whence 0 = ^(1“- = 1 - f m*^, 

Z¥ 


A = 
5 = 


4(1-4) 

- Sm" 
1-4 


= - 

= OT"*, 


^ - 2 {1 - (1 - 2m + f m^} “ 


D = 


3mV 


= - f 


2 (1 — m'*) 

^- 4Tl-l-2m + ^ 

Therefore 

/I — 

+ e cos {cd — a) 

— cos 2 ( 90 -~y) 

u ~ ai^ + cos {(2 — 2m) 6 — 2^} ^ 

+ ^me cos{(2-2??i--c) 0 - 2/3 + a} 

— |mV cos(m0 + yS — ?) 

1^ + cos{(2 - 2?w - 2c) 0 - 2/9 + 2a} J 

49 The differential equation of the latitude 

d^s Ps-- 8 T ds /d^s \ r T 
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becomes, aftei substitution, 


(fs 

W 


+ 5 = 


Assume s — 


— fmVc sm {gd — 

-I- sm {(2 ~ 2m — — 2/3 H- 7} 

h sin {gO — 7) 

+ A sm{(2 — 2m — 0 - 2/3 + 7}- 


Then, by substitution, we get 

7c (1-/) =-|m^^, 

A {1 - (2 - 2m - gY] = + 
Theiefoie = V(H“ = 1 + 


therefore 



2 {1 - (1 — 2m — ) } 

A:sm(5'^ — 7) 

+ §mA: sm {(2 — 2m — 0 — 2^ + 7} 


4. 


60 We can now find the connexion between the longitude 
and the tune to the second order, 



and from Art (48) we have 

1 1 

hu^ haf {1 + e cos(c0 — a) + + S,}" ’ 

hemg the sum of all the terms of the second, and 
those of the third order m u,) 


= i [1 - 2 (e cos[c9 - a) + + S^} 

+ 3 {e cos(c0 — a) + Sg + Sj* — &c.] 

= [1 — 2e cos {c0 — a) — 2Sj, — 22, 

** + 3e“ cos* {c9 - a) + 6S,e cos {cB- a)] 
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— ^ [1 — 2e cos {od — a) — 2S^ — 2S3 + cos2 {c0 — a) 

4 cos {(2 — 2m — 2c) 0 — 2jS 4 2a}] 

1 4 Ic' 4 f 4 m"* — 2e cos [c6 — a) 

4 fe” cos2(c0 — a) 

4 cos2(^0 — 7) 

— 2m“ cos {(2 — 2m) 0 - 2^8} 

— i^mccos{(2— 2m-c)8-2/34a} 
4 3mV cos (m8 4 y8 — 

— Vme'*cos{(2— 2m— 2 c) 8— 2;842a} 

Also, fiom Art (45), 

r T 

1 — \j^ dd = l - |m® cos {(2 — 2»j) 0 — 2/3} 

4 cos{(2-2m— 2c) 8 - 2/J 4 2a}, 

neglecting the other term of the third ordei, the coefficient of 
the argument not bemg small 

We have now to multiply these results together, and we see 
that the term havmg for argument (2 - 2m - 2c) 8 - 2^ 4 2a 
will disappear in the product If we tiace this term, we shall 

r T 

find that it aiose in from letaming originally terms of 

the fourth order, but in it aiises from combining terms 

oiigmally of the fiist and third orders If, thcrefoie, we had 
1 ejected terms beyond the thud order indiscriminately, the ex- 

piession ^ would have contained this term, intioduced by ~ , 

to the third order, and m t it would have been raised to the 
second order, and theiefore formed an impoitant part of its 
value instead of disappearmg altogether from the expression 
Hence the necessity for retaining such tenns of the foiiith ordei 
T 

1 “ 12 -# 

« M 
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1 + + vi^ — 2e cos (60 — a) 

+ \e‘ cos2(c0— a) 

^ cos2 {gO - 7) 

di) ha^ \ — cos{(2 — 27?^) 6 — 2/3} 

— ^-^me cos {(2 — 2m — c) 0 — 2^ + a} 
+ cos [mO + /3 — f) 

Let + + 

therefore = 1 - 4e" - 4^^ - m' to the third order , 

ha 

therefoie, multiplying by p and integrating, we get, still to the 
second oidei, 

jpt=0'-2e sin {cO — a) + sin2 (c0 — a) 

+ \'k^ sin2 (9^0—7) 

— y m‘^ sin {(2 — 2m) 0 — 2/3} \ 0^, 

— sin {(2 — 2m — c) 0 - 2/3 + a} 

+ 3me sm (m0 H- iS - ^) 

no constant is added, the time being reckoned fiom the instant 
when the mean value of 0 vanishes, for the reasons explained 
in Art (34) 

61 The pieceding equations C/,, 0^ give the leciprocal 

of the radius vector, the latitude and the time in terms of the 
tiue longitude, but the principal object of the analytical Inves- 
tigations of the Lunar Theory being the formation of tables 
which give the coordmates of the moon at stated times, we must 
express 5, and 0 in terms of t. 

To do this, we must reverse the senes = 0 — &c , and 
then substitute the value of 0 in the expressions for u and s. 

0 =rp + 2e sin(c0-a) to the fiist order 

— pt 2e sin {cpt---a) . • • 5 


Now 



46 


LUNAR THEORY 


therefore cO — a cpt — a + 2 e sin(cp^ — a) to the first order, 

2e sm(c6 - a) = 2 e {sm(cp< — a) + 2 e sm (cp^ — a) cos{q)t- a)} 

to the second order, 
= 2 e sm[cpt— a) + 2e® sin 2 (cpt— a) , 

and as 9 and pt diiSPer by a quantity of the first order, they may 
be used indiscriminately in terms of the second order , therefore 

9 —pt + 2e B>m{cpt — a) + f sm 2 [cpt — a) 

— sm2 {gpt — 7) 

+ y sm{(2 — 2w) pt — 2 ^] 0'^ 

+ ymesin{(2 — 2 w2-“c)^^— 2)S -f a} 

— 3W sm [mpt + — S') . 

52 In the value of u given in Art (48), substitute pt for 
9 in terms of the second order, and pt 4 2 e sin (cp^ — a) m the 
term of the first order , then 

1 — |F — — + e cos(cp^ — a) + ^ cos2 {cpt — cf) 

-\k‘ cob 2 {gpt - i) 

+ cos {(2 - 2 m) pt — 2/3} 

+ ymecos((2— 2m— c)j9^"-2/3+a} 

the other terms m the value of u m Ait (48), which were there 
retained only foi the sake of subsequently finding being of 
the third ordei, aie here omitted 

53 Similarly, the expression for 5 becomes 
s = L sin{(p_pi — 7) 4* 2 e sin(cp^— a)} 

4 %m 1 c sm {(2 — 2 m — 5^) — 2^ 4 7) , 

sm(9p^ — 7) 

4 e sm [[g 4 0) — a — 7} 

— e sm {(^r — c) _p« 4 a — 7} 

4 fw sin{(2 — 2m-9)p?^ - 2/3 4 7} 

The expression for s is more complex in this form tluui when 
given in terms of the true longitude 9 
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54 If P be the moon’s mean paiallax, and n the parallax 
at the time 

" = S^ceof 2) = r-, : = (1 - third order, 

-V(i + 0 

= Ru {1 — cos 2 {gjpt- 7)} to the second order, 

( 1 — 4- e coB[cjpt — a) + e** cos 2 [cpt — a) 

+ w® cos {(2 — 2m) — 2yS} 

+ cos{(2— 2?}^ — c)j9^ — 2/3 + a}, 

but P = the portion which is independent of periodical teims, 

= Ea (1 — Ti? — \in^ - e®) ; 

( 1 4 e cos(cp^ — a) + cos2(cpi — a) 

4 m** cos{(2 — 2m) — 2y8} 

4 cos{(2 — 2m — c) — 2/8 4 a} 
neglecting terms of the third order 


55 Here we tenninate our approximations to the yalues 
of w, 5, and 0 If we wished to cany them to the third order, 
it would be necessaiy to include some terms of the fourth and 
jfifth orders according to Art (29), and the values of P, P, 
and 8, given m Art (23), would no longer be sufficiently 
accurate, but we should have to recur to moie exact values, 
and from them obtain terms of an order beyond those abeady 
employed 

If this be done, it is found that 


P 


a n + 


_ R — M „a 


}. 


T 

h‘u^ 


E-M 

jEJ 4- Jf 


— sin (0 — 0') 
a ^ ^ 


These terms of the fourth order become of the third order in 
the value of w, and therefoie also of t, the coefficient of 0 being 
near unity 

We shall see further on (Appendix, Art 97), to what purpose 
a knowledge of the existence of these terms has been apphed. 
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56 The process followed m the pieceding pages is a suffi- 
cient clue to what must be done foi a higher appi oximation 

The coordinates vl and ff of the sun’s position are, by the 
theory of elliptic motion, known in terms of the time and t 
IS given in terms of the longitude 6 by the equation Hence 
v! and & can be obtained in teims of 0 , but it will be necessary 
to take into account the slow progressive motion of the sun’s 
perigee, which we have hitherto neglected This will be done 
by writing c' 6 ' - ^ foi & - f, d being a quantity which diffeis 
very little from unity ^ 

These values of 0 ', together with those of u and s m terms 
of 0 , as given by and 8^^ are then to be substituted in the 
coiTected values of the forces, and thence in the differential equa- 
tions The mtegiations being performed as before will give the 
values of w, 5 , and t m teims of 0 to the third older, and from 
these, as m Aits (51), (52), and (53), may be obtained w, 5 , and 0 
in teims of t 

57 More appioximate values of c and g are obtained at 
the same time, by means of the coefficients of cos (c 0 — a) and 
sin (^ 9 ^ 0 - 7 ) in the differential equations, (see Appendix, Aits 
94 and 95) 


* ‘ En reflSchissant sur les termes quo doivent introduire toutes les quan- 

* tit6s prec6deiites, on voit qu’ il se pent glisser des cosmus de Tangle 6 dont 
‘ nous avons vu le dangereux effet d'amener dans la valeui de u des arcs au 
*lieu de leuis cosmus, de tels termes viendiont, par exemple, de la com- 

* binaison des cosmus de (1— w)6 avec des cosmus de m0 

* Pour e-vuter cet inconvenient qui Oterait k la solution preeddente 

* Tavantage de convenir aun aussi grand nombre de revolutions qu’on voudrait, 
‘ et la pnverait de la simplicite et de Tuniveisalite si precieuses en matlioma- 

* tiques, il faut commencer par on chercher la cause Or, on decouvre facile- 
‘ meiit que ces termes ne -viennent que de ce qu’on a suppose fixe Tapogee 
‘ du soleil, ce qui n'est pas permis en toute ngucur, puisquo quelquo petite 
‘ que soit sur cet astre Taction de la lune, elle n’en cst pas moms recllc ct 

* doit lui produiie un mouvement d’ apogee quoique tres lent u la veritc ’ 

Clairaut, TheoriQ de la lAine^ p 55 2me Edition 
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58 The values to the fourth order are then obtained fiom 
those of the thud by continuing the same process, and so 
on to the fifth and higher orders, but the calculations are so 
complex that the approximations have not been carried beyond 
the fifth order, and already the value of 0 in terms of t contams 
128 periodical terms, without including those due to the dis- 
turbances produced by the planets The coefficients of these 

periodical terms are functions of e\ “ , c, A, and are 

Cb 

themselves very complicated under their hteral forms, that of 
the term whose aigument is twice the difference of the longitude 
of the sun and moon, for instance, is itself composed of 46 terms, 
combmations of the preceding constants. 

See Pont^coulant, Systlme du Monde^ tom iv p. 572, 
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CHAPTER V 


NlIMEEtCAL VALUES OP THE COEFFTCIEKTS 

ijO Having thus, fiom thcoiy, obtained the fimii of the 
devolopincnts of the cooidiuates of the moon’s position at any 
time, the next neccssaiy stop is the deteimination of the 
iiumciical values ot the coefficients of the seveial terms 

We lieic give three ditfeicnt methods which may ho em- 
ployed foi that piupose, and those may, moieover, be combined 
act 01 ding to ciicvnnstanccs 

()() fi'iiht mcfliotJ By paiticulai obseivations of tlic sun 
and moon (/ c by obseivations made when they occupy parti- 
(ular and selected positions), and also by obseivations scpaiatcd 
by vi'iy long intoivals, such, for instance, as ancient and modem 
eclipses, the values ot the constants jo, wi, a, / 3 , 7, f, which enter 
into the and of the additional ones which enter into 

the coefficients of the teims in the picvious dovclopincnts, may 
Ik> obt.iined witli great imeiiiacy, and by thou means, the coeffi- 
cients themselves; c and ff being also known in terms of the 
othoi constants 

I’liese may piopeily be called the theoretical values of the 
i-oeffiiients, the only lecouibo to obseiwation being for the detei- 
nnuation ot the nunieiical values of the elements 

(51 Hrcoml ‘method Let the constants which enter into the 
aKIiimeuh bo deteimincd as in the fiist method ; and lot a large . 
number of obseivations bo made, fiom each of which a value 
of the tiuc longitude, latitude, or parallax is obtained, together 
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with the coi responding value of t leckoned from the fixed epoch 
when the mean longitude is zero Let these corresponding 
values be substituted in the equations, each observation thus 
giving use to a relation between the unknown constant coeffi- 
cients. 

A very great number of equations being thus obtained, they 
are then, by the method of least squares or some analogous 
process, reduced to as many as there are coefficients to be 
detci mined The solution of these simple equations will give 
the requiicd values 

This method, however, would scarcely be practicable in a 
high oidci of appioximation For instance, in the fiftli order, 
as stated m Art (58), each of the numerous equations would 
consist of 130 terms, and these would have to be reduced to 
129 equations of 130 terms each 

62. Third method When the constants which enter into 
the aiguments have been determined by the first method, we 
may obtain any one of the coefficients mdependently of all the 
others by the following piocess, provided the numbei of obseiv- 
ations be very great 

Let the foim of the function be 

F == ^ -1- jB sm^ + G smef} h- &c , 

and let it be required to determine the constants A, j?. O', &c. 
sepal atcly ; 0, &c being functions of the time. 

Lot the icsults of a great number of observations corre- 
spondmg to values d^, &c , (f>^, ^ 3 , &c , be F„ F,, F„ 

«&c ; so that 

V^ = A -{■ S sin^i + G am<f>^ + &c , 

V^=A-\-B sm ^3 + 0 + &o., 

V^ — A + B + C sm ^3 + &c , 

A + B sin0„ + G + &c. 

Now, n being very great, we may assume that the periodical 
terms will, on the whole, be as often positive as negative ; and 

£ 2 
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tliereforej that if we add all the equations togethei these terms 
will cancel one another ; 


therefoie 


^ ^ ^ 


+ f: 


which determines the non-periodic part of the function 

To determine B Let the observations be divided into two 
sets separating the positive and negative values of sin0; then 
the other periodical terms, not havmg the same period, may be 
considered as cancelling themselves in adding up the terms of 
each set Let there be r terms m the first set and s terms m 
the second, and let F', F*" be the values of F corre- 

sponding to positive values of sin0, which values we may assume 
to be uniformly distributed from sinO to sinw, and therefore 
to be smSdj sm2S0, sinr 80, where r 80 = tt 

And, again, lot F], F],, F]^,, , F^, be the values of F corre- 

sponding to the negative values of sm0, viz . sin(— A0), sm(— 2A0), 
sin(— s A0), where s A0 = tt Then, 


F'=-4 + 5sin80 +Osm«^>'H- , 

F" = ^ + i?sin2.S0 + (7sin<^"-i- 


V'" =:A + B sinr Bd +0 + ; 

therefore 


— J3sin5.A0 4-(7sm<)b, -h ; 
therefore 

F, + F,+. +V,^sA-B2;{Bme) 


= rA+^ HJi&md.Be)-, 


therefore 


»• '^j» 


d0 


~ A « 

TT 


F^ — j5sinA0 -f-C7sin<^^+ ,, 
F^, = -4 — 5sm2 A0H-(7sin^^^-f , 




B 


therefore 




Z±L± ±L^A-^f,^ede 

2B 


= A-- 


TT 


therefore B = 


TT /V’+V"+.. V' 
4 I 




T 


a 
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and in a similar maimer may each of the coefficients be inde- 
pendently determmed * 

For further remarks on this method, see Appendix, Art (104) 


* If r and 8 are not sufficiently great to allow us to substitute JJ' smOdO 
for sm0 80, we must proceed as follows 

V' + V" + -f F*" = r-4 H- ^ (sin 50 + sm250 + + sinr50) 

. ^ sin H- 1 ) 80 sin |r 80 

= rA + M — ^ ^ 

8mi86 


= rA + B 


sm i50 * 


r + nearly 

r r iS9- /gae* ’ 

1 -J- - 




Similarly, 
therefore B s 




_ 1 


8 


la ^j, 


ir 

/r'+F"4 

+r' r,+r„+ 

+^A 

■‘{‘-hI 


8 


H‘*£' 

/I , iw /F'+r'+ 

\r* 5*/ J \ > 

+ F»' 

8 

4F,^ 
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CHAPTER VI. 


PHYSICAL INTERPRETATION 

63 The solution of the problem which is the object of the 
Lunai Theory may now be considcicd as effected; that is, we 
have obtained equations winch enable us to assign the moon^s 
position m the heavens at any given time to the second order 
of approximation ; wo have explained how the numerical values 
of the coofHcicnts in these equations may be dctei’mincd fiom 
obseiwation ; and we have, moicovci, shewn how to proceed m 
Older to obtain a highoi appioxiination ^ 

It will, however, be interesting to discuss the icsults we have 
arrived at, to see whether they will enable us to fom some 
idea of the nature of the moon’s complex motion, and also 
whether they will explain those inequalities oi departuies fiom 
lunfonn ciiculai motion which ancient astionomcis had obseivcd, 
but winch, until the tunc of Newton, were so many unconnected 
phonoincna, or, at least, had only such arbitrary connexions as 
the astionomers chose to assign, by giafting one eccentnc or 
epicycle on another as each newly discoveicd inequality seemed 
to render it necessary. 

It is true that our expressions, composed of periodic terms, 
arc nothing more than tianslations into analytical language of 
the epicycles of the ancients ;t but they aie evolved directly 


* The mcfULs of talcing into account the ellipsoidal figure of the earth and 
the disturbances produced by the planets, axe too complex to toi-m part of an 
introductory treatise Por information on these points reference may be made 
to Airy's Fiqiiro aftko Uarih Ponte coul ant's Si/stdme du MondCy* vol iv 
t Sec Whcwcll's “Uistojy of the Inductive Sciences " 
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from the fundamental laws of force and motion, and as many 
new terms as we please may be obtained by carrying on the 
same process; whereas the epicycles of Hipparchus and his 
followers were the result of numerous and laborious observa- 
tions and comparisons of observations , each epicycle bemg 
introduced to correct its predecessor when this one was found 
madequate to give the position of the body at all times: just 
as with us, the terms of the second order correct the rough 
results given by those of the first, the terms of the thud order 
correct those of the second, and so on But it is impossible to 
conceive that observation alone could have detected all those 
minute irregularities which theory makes known to us m the 
terms of the third and higher orders, even supposing our in- 
struments far more perfect than they are, and it will always 
be a subject of admiration and surprise, that Tycho, Kepler, and 
then predecessois should have been able to fad their way so far 
among the Lunar inequalities, with the means of observation 
they possessed 


LONGITUDE OF THE MOON 

64 We shall fiistly discuss the expression foi the moon’s 
longitude, as found Art (51) 

0 = + 20 sm [cjpt — a) -f sin 2 — a) 

+ sin {(2 — 29W — c) — 2y8 H- a} 

4* sm{(2 -2^)^^- 2/9} 

— Zmd sin + /3 — S’) 

— Bm2(;f7j?« — 7) 

The mean value of 0 is ; and m order to judge of the effect 
of any of the small terms, we may consider them one at a time 
as a correction on this moan value or we may select a com- 
bination of two or more to form this correction* 

Wo shall have instances of combinations m explaining the 
elliptic inequality and the evection^ Arts (66) and (70) ; but m the 
remaining inequalities each term of tlie expression will form a 
conection to be consideied by itself 
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65 Neglecting all the periodical terms, we have 

ae 

which indicates uniform angular velocity; and as, to the same 
order, the value of u is constant, the two together indicate that 
the moon moves uniformly m a circle, the period of a revolution 
Stt 

being — 5 which is, therefore, the expression for a mean sidereal 

month, or about 27-J- days.'* 

The value of^ is, according to Art (50), given by 

and as m is due to the disturbing action of the sun, we see 
that the mean angular velocity is less, and therefore the mean 
periodic time greater than if there were no distui banco 


Elliptic inequality ot Equation of the Centre, 

66, We shall next consider the effect of the first three terns 
together* the effect of the second alone, as a coriection of 
will be discussed in the Histoncal Chapter, Art, (109) 

0 + 2e sm(c^i- a) + sin2 a), 

which may be written 

0 + 2e sin [pt - {a + (1 - c) pt]] + sm2 |>j5 - {a + (1 -c)pt]'\. 

But the connexion between the longitude and the time in an 
ellipse described about a centre of foicc in the focus, is. Art, (13), 
to the second older of small quantities , 

6 = nt-\-2e sm(n^ — a') + sin 2 (?^^ — a'), 


* The accuiatc value was 27d 71i ISm 11 261s m the year 1801 See 
Art (99) 
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where n is the mean motion, e the eccentncity, and a' the 
longitude of the apse.* 

Hence, the terms we are now considering mdicate motion 
in an ellipse ; the mean motion hemg the eccentricity e, and 
the longitude of the apse a + that is, the apse has a 

piogressive motion in longitude, uniform, and equal to (1 - c)^ 

67. The two terms sin(c^«- a) + sin2(c^f-a) con- 
stitute the dlvptic, inequality^ and their effect may be farther 
illustrated by means of a diagram 

Lot the full line A.MB (fig 8) represent the moon’s orbit 
about the earth JE^ when the time t commences, that is, when 
the moon’s mean place is m the prime radius -E^T from which 
the longitudes are reckoned 

The angle TEA^ the longitude of the apse, is then a At 
the time when the moon’s mean longitude is TjEM* = the 
apse line will have moved in the same direction through the 
Jingle AEA' = (] - c) TJEM, and the orbit will have taken the 
position indicated by the dotted ellipse , and the true place of 
the moon in this orbit, so far as these two terms aie concerned, 
will be where 

MEm = 2e sm [c^t — a) + fe** sin 2 [cjpt — a) 

= 2e suiA'EM + mi2AEM 
:r: 2e mxAEM{l + |e c.o^AEM) , 

which, smce e is about is positive from perigee to apogee, 
and therefore the true place before the mean^ and the contrary 
from apogee to perigee : at the apses the places will comcide 

68. The angular velocity of the apse is (1 — c)y>, or, if for 
c we put the value found m Art (48), the velocity will be 
Hence, while the moon describes 360'’, the apse should describe 

360“ = nearly, m bemg about ^ 


* The epoch e which appears m the expression of Art (13) is here omitted , 
a proper assumption for die origin of t, as explamedm Art (34), enablmg us 

to avoid the e 
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But Hipparchus had found, and all modem observations 
confirm his result, that the motion of the apse is about 3° in 
each revolution of the moon See Art (112) 

This difference arises from our value of c not being repre- 
sented with sufficient accuracy by 1 - 

Newton himself was aware of this apparent disci epancy 
between his*tlieory and observation, and we are led, by his 
own expressions (Scholium to Prop 35, lib iii in the first 
edition of the Princi^ia)^ to conclude that he had got over the 
difficulty This is rendered highly probable when we con- 
sider that he had solved a somewhat similar problem m the 
case of the node; but he has nowheic given a statement of 
his method and Clairaut, to whom we aie indebted for the 
solution, was on the point of publishing a new hypothesis of 
the laws of atti action, in order to accoimt lor it, when it oc- 
curred to him to caiTy the approximations to the thud ordei, 
and he found tlie next term in the value of c nearly as con- 
sideiable as the one already obtamed See Appendix, Art (94). 

c = 1 — j 

1 — c = (1 -J- ’'^m) , 

• (1 - c) 360° = (1+ /xrj) (value found previously) 

= 2f° nearly, 

thus reconciling theory and obseiwation, and removing what had 
proved a great stumbling-block m the way of all astronomers * 
When the value of c is earned to higher ordois of ap- 
proximation, the most perfect agreement is obtamed 

The motion of the apse line is considered by Newton m his 
Prmcijpzaj lib I , Prop 66, Cor. 7. 


Pvcction 

* 

69. The next term in the value of 0 is 

H- ^me sin{(2 — 2m — c)^^- 2/3-l-a)} 


* See Dr WliewcH'R B) idf/eivato 
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We shall considei this term in two ditFerent ways. 

Firstly, by itself, as forming a correction on 

d ^ jpt + sin{(2 — 2?^ — — 2yS-l- a} 

Let 3) = ]^t = moon’s mean longitude at time 

O = mjpt + y8 = sun’s . . , 

a' = (1 — -H a = mean longitude of apse , 

then 

0 = sm[2{j9^— ^)] — [pt— (1 — c)pt-{' a]] 

— pt + sm[2 ( }) “ o) — ( 2) - a')] 

The effect of this term will therefore be as follows . 

In syzygies 0 — '^me — a') ; 

01 the true place of the moon will be before or behind the 
mean, according as the moon, at the same time, is between 
apogee and perigee or between perigee and apogee. 

In quadratures Q = pt ^ ^me sin ( J - a'), 

and the circumstances will be exactly reversed 

In both cases, the correction will vanish when the apse 
happens to be in syzygy or quadrature at the same time as 
the moon 

In intermediate positions, the nature of the correction is 
more complex, but it will always vanish when the sun’s lon- 
gitude IS equal to the arithmetical mean between those of the 
moon and apse, or when it differs from it by any multiple 

of 90'* ; for if O = — 'r 90°, 

sin[2 ( 5 — O) — ( J) - a')] = sm ( 3) -1- a' — 2o) 

— sinr 180° 

= 0 . 

70. The other and more usual method of considering the 
effect of this term is in combmation with the two terms of 
the elliptic inequality, as follows: 

To determine the change in the position of the apse and in 
the eccentricity of the moon^s orhit pi oduced hy the ejection. 
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Taking the elliptic inequality and the evcction together, 
we have 

0 = jpt -h 2e sm [cpt — a) + sin 2 {cpt — 

+ ^me sin {(2 — 2m - c)pt — 2;8 4* a} 

Let a! be the longitude of the apse at time t on supposition 
of uniform progression, 

O . . . ... sun : 


whence a' = (1 — c)pt + a, 

O =!= mpt 4- (3 
-Ajud the above may be written 

^ 4 2e sin [cpt — a) 4 sin 2 [(pt — a) 

4 ^me sm [cpt - a 4 2 (a — O) J , 
«ind the second and fourth terms may be combined into one, 
2E sin(cp^ — a4S), 

if -B cosS = c 4 cos 2 (a' - O), 

E sinS = sin 2 (a' — O) , 

whence tan S = ^ 

1 4 cos 2 (a — o) ^ 

E^ = 4 cos 2 (a' — O)}® 4 <i^[ sm 2 (a — O)]'*^ ; 

or, approximately, 


S = ^^7)1 sin 2 (a' - O), 

^ = ^{ 1 - 1 - x^ni cos2(a' ~ o)). 

The term sm2(cpi — a) will, thoioforc, to the hccoiuI 

Older, be expressed by 

1 ^-* sm2(cpi5-a4 S), 
and the longitude becomes 


6 =:pt + 2Esm{cpt - a 4 S) 4 l;E^ sm 2 {(pt - a 4 S), 
or d==pt+ 2E&m{pt-(x:+S) 4 sm2(;^^- a'4 S) ; 

bnt the last two terms constitute elliptic inequality m an oibit 
whose eccentricity is E and longitude of the apse a' - S ; there- 
fore the evection, taken in conjunction with elliptic niequahty, 
has the effect of lendering the ccccntucity of the moon’s orbit 
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vanable, increasing it by when the apse-line is in syzygy, 
and diminishing it by tie same quantity when the apse-lme is 
in quadrature ; the general expression for the increment being 

^me cos 2 (a' — O) 

And another effect of this term is, to diminish the longitude 
of the apse, calculated on the supposition of its uniform progres- 
sion, by the* quantity 8 = sin2(a — O) ; so that the apse 
IS behind its mean place from syzygy to quadrature, and before 
it from quadrature to syzygy* 

The cycle of tfeese changes will evidently be completed m 
the period of half a revolution of the sun with respect to the 
apse, or m about of a year 

,71. The period of the evection itself, considered mdepen- 
dently of its effect on the orbit, is the time in which the 
argument (2 — 2m — o) — 2/3 + a will increase by 27r 
Therefore penod of evection 

27 r mean sidereal month 


(2 — 2m — 2 — 2m — c 



= 31-J- days, nearly .f 


Newton has considered the evection, so far as it arises from 
the central disturbmg force, in Prop. 66, Cor, 3 of the PmVz- 
cvpia 

Variation 

72 To exjplain the physical meaning of the tenrm 
y m'* sm {(2 — 2m) pt — 2)8}, 

%n the eapression for the moorCs longitude^ 

0^pt^ ym'® sm{(2 — 2m)j?^ — 2)8}. 


* The change of eccentncity and the yanation in the motion of the apse 
follow the same law as the abscissa and ordinate of an ellipse referred to 
its centre for if — a = a; and 5 = y, then 

V* 

t The accurate value is 31 8119 days 
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Let }) represent the moon’s mean longitude at time 
O sun’s 

therefore D 

O = 7}ijpt -I- /3 , 

and the value of 6 becomes 

^=^9^4- ym^sin2( 3) — o), 

which shews that fiom syzygy to quadiatme, the" moon’s true 
place is before the mean, and behind it fiom quadratuie to 
syzygy; the maximum diffeicnce being in the octants 

The angular velocity of the moon, so tar as this term is 
concerned, is 

dO 

^ — cos2(3) — o), 

—p {1 + cos 2 ( 3) ~ o)}, nearly, 
which exceeds the mean angular velocity p at syzygies, is equal 
to it m the octants, and less in the quadiatuies 

This inequality has been called the Variation^ its peiiod is the 
time in which the argument (2 — 2 w 2 ) pt — 2y8 will increase by 27 r , 

penod of variation = , = mean synodical month 

^ (2 — 2?w) p 2 

’ = 14f days, nearly.* 

73 The quantity is only the fiist term of an endless 
series which constitutes the cocflScient of the vaiiation, the other 
terms being obtained by caxrymg the approximation to a higher 
order. It is then found that the next teim m the coefficient 
IS which IS about of the first term ; and as there are 

several other important terms, it is only by carrying the ap- 
proximation to a high order (the 5th at least) that the value of 
this coefficient can be obtained with sufficient accuracy from 
theory In fact, ig’-m’* would give a coefficient of 28' 32" only, 
whereas the accurate value is found to be 39' 30". 

The same remark apphes also to the coefficients of aU the 
other terms 


The accurate value is 14 766294 days 
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74: As fai as teims of the second order, the coefficioi^f 
the variation is independent of e the eccentricity, and. tlio 
inclination of the orbit It would therefore be the sstxxi® 
an orbit onginally circulai, whose plane coincided with the 
plane of the ecliptic it is thus that Newton has considoirod it- 
Pnnc'ijf, Prop 66, Cor 3, 4, and 5 


Annual Equation 

75 To explain the physical meaning of the term 
— 3W sin [mpt-^ 

in the expression for the moorCn longitude 
6 = pt - ^m{mpt + /9 S’)? 

= ^^-- 3W sm (longitude of sun - longitude of sun’s pex-igec), 
^ jpt Zmd sin (sun’s anomaly) 

Hence, while the sun moves from perigee to apogee, the true 
place of the moon will be behmd the mean, and from a^pogee 
to perigee, before it The period bemg an anomalistic year, 
the effect is called Annual Equation, 

Differentiatmg 9 we get 
dO 

~ = ^[1 ~ 3m V cos (sun’s anomaly)} 

Hence, so far as this inequality is concerned, the moon’s an- 
gular velocity IS least when the sun is in perigee, that is eit 
present about the 1st of January, and greatest when tlxo stiii 
IS in apogee, or about the 1st of July 

The annual equation is, to this order, mdependent of the 
eccentricity and inclmation of the moon’s orbit, and tlioircforo, 
like the variation, would be the same m an orbit orig-inally 
cncnlai Vide Newton, Principia^ Prop, 66, Cor. 6 

Reduction 

76. Before considering the effect of the term — — 

which, as we shall see Art. (82), is very nearly equal -to tlie 
difference between the longitude in the orbit and the longitude 
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U1 the ecliptic, it will be convenient to examine the expression 
or t e atitude of the moon, and to see how the motion of the 
node IS connected with the value of 

LATITUDE OF THE MOON 

77 . The expiession Art ( 49 ), found for the tangent of the 
latitude,* IS 

s = h sm(5'0 - 7) + Imk sm{(2 - 2w — 0 - 2/3 + 7} 

If we reject all small terms, we have 

s = 0, 

or the orbit of the moon comcidmg with the echptic, which is 
a first rough approximation to its time position 

78 Takmg the fiist term of the expansion 

9 = Te sm(^0 — 7) 

we may write it 

s = ^8m[0-{7-(_^_ 1) 0}] 

Let 'TNm be the ecliptic (fig 9), N the moon’s node when 
her true longitude is zero, and let M be the position of the 
moon at time «, m her place referred to the ecliptic; 
therefore Ti7=7, Tm = 0 , Un Mm = s 

Take NM = (,9 - 1) 0 m a letrograde direction, and jom MJST by 
an arc of great circle , 

then sin N'm = cotilfiVW, 

or sm [0 - {7 - (5- _ 1) 0}] = s cx>iMN'm , 

which, compared with the value of $ given above, shews that 
MN'm = tan"’^ is constant, and therefore the term k sm{ffd — 7) 
indicates that the moon moves m an orbit inclmed at an angle 
tan''^! to the ecliptic, and whose node regredes along the ecliptic 

with the velocity (5' — 1) ^ > or 'with a mean velocity (9 — 1) 


• This expression for the tangent of the latitude is more convenient 
that which gives it in terms of the mean longitude, Art (53) on account of 
tiie loss niinibor of terms involvedt See Pont^coulent, vol ly j p, 630* 
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79. 


Hence the 


period of a revolution of the nodes = 
_ one sidereal month 


2‘7r 

{9 - l)i’ 


9-1 

but, from Art (49), the value of 9 = 1 + fm* ; 
therefore penod of revolution of nodes = sidereal month 


= 6511 days, nearly 


This will, for the same reason as in the case of the apse, 
Ai*t (68), be modified when we carry the approximation to a 
higher degiec, tins value of g is, however, much more accurate 
tlian the coiiesponding value of c, for the third term of g is 
small; the value to the third order hemg (see Appendix, 
Alt. 95) 


= 1 + I m" - 


and the period of revolution of the ^odes = 


one sidereal month 


fm^(l — fm) 

= 6705 days, neaily- 


This IS not far from the accurate value as given by obser- 
vation, and when the approximation to the value of g is carried 
to a higher ordci, the agieement is nearly perfect 

The true value is 6793 39 days, that is about 18 yrs , 7 mo 


JEvection %n Latitude 

80. To explain tJie variation of the inclination and the irre^ 
gularity in the motion of the node expressed hy the tefrm 
+ %mh sin {(2 — 2m— g) 0 — 2/3 + 7}. 

This term, as a correction on the preceding, is analogous to 
the ovection as a coirection on the elhptic mequality. 

Taking the two terms together, 

s mi[gd — 7 ) + %mTc sin {(2 — 2 m -^) 5 — 2 ^ + 7 } 

Lot D = longitude of moon = 

0 =.***. •• sun =^^ + /S, 

12 = . . node = 7 — [g — \)6^ 

s sin(2) — i2) •\-%mh sm{3)-X2 — 2(0— X2)} 

F 


therofoie 
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Now these two teims may be conibmcd into one 
s = K sin( ]) — 12— S), 
if ^cosS = i cos2 (o — 

K sin S = I mli. sm 2 ( O — 12) , 

siii2(o — 12) 


whence tan S = 


1 + I'iii cos2(o — 12) ’ 


= A.^{1 + >cos2(o-i2)r + >’'11 2(0- 


oi approximately, 

S = sin 2(0 — 12), 

K = cos 2(0 — Q<) } , 

but the equation 

s — K sin( 3) —12— S) 

lepresents motion in an oibit inclined at an angle tan^Vr to the 
ecliptic, and the longitude of whose node is 12 + S. 

This teim has thciefoic the following effects 
1st. The inclination of the moon’s oibit is vaiiablo, its tan- 
gent increases by when the nodes aic in syzygics, and 
decreases by the same quantity when they aie m quadrature; 
the geneial expression for the incieasc being cos 2(o— 12). 

2nd The longitude of the node, calculated on supposition 
of a uniform regression, is incicased by S = sm 2(0—12), 
so that the node is bcfoie its mean place while moving fiom 
syzygy to quadrature and behind it tiom quadratuie to syssygj^ 

PHncii)ia^ book III , props 33 and 35. 
The cycle of these changes will be completed in the poiiod 
of half a revolution of the sun with icspect to the node, that is, 
in 173 21 days, not quite half-a-year 


81 The tangent of the latitude has heie been obtained ; if 
we wish to have the latitude itself it will be given by the 
formula 

latitude = s — Js® + -^ 5 ® — &c , 

which, to the degree of approximation adopted, will cloaily bo 
the same as s. 
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Reduction 


82 We may now consider the teim which we had neglected 
Art 76) m the expiession foi the longitude, namely, 

— ^ 1 “ sin2(^j)i— 7). 

Let A7 (fig 10 ) be the position of the node when the moon’s 
ongitude is 0, M the place of the moon, m the place referred 
0 the ecliptic. 

Therefore Tm = 5 , 

TA7=7- (^-1)^, 

Nm gQ — 

imN = h 


The nght-angled spherical triangle NMm gives 


cosAT = 


tanNm ^ 
tanATi/’ 


herefore 1 — cos A^ tanA/'if — tmNm 

1 + cos A^ tanA^if 4 - twnNm ’ 

'T tan”— = 

2 sm(A7M'+ A7«i) ’ 

r, smce both N and NM- Nn arc sm all, 

tanW NM-Nm 
~ “ ~1^2Nm aPPros^ately ; 

herefore NM-Nm = sm 2 (.^ 0 - 7 ) = sin 2 (.^ 2 i<- 7 ), nearly 

Hence this term, which is called the reduction^ is appioxi- 
aately the difference between the longitude in the orbit and 
he longitude in the ecliptic 


RADIUS VECTOR 

83 To explain the physical meaning of the terms in the 
aim of u. 

We shall, for the explanation, make use of the formula 
rhieh gives the value of « m terms of the true longitude, 
urt (48). ^ ’ 
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Firstly, neglecting tlie periodical terms, -we have for tlic 
mean value 

The term - \m% which is a consequence of the distuihing effect 
of the sun, shews that the mean value of the moon’s radius 
vector, and therefore the orbit itself, is largei than if tlicie 
weie no disturbance 


JSlhptie Inequality 

84 To explain the effect of the term of the jSnxf ordor^ 
u = a{l + ficos(c0 — a)}, 

= a[H-ecos0 — {a+ (1 — c) 0]] 

This 13 the elliptic twguahty, and indicates motion ui an 
ellipse whose eccentiicity is e and longitude of the apse 
a 4- (1 _ c) 0 , and the same conclusion is drawn witli respect 
to the motion of the apse as m Art (66) 


Evectio7i 

85 To ex^lmn the jphysical meaning of the term 
^mea cos{(2 - 2in — c) 6 - 2^ 4- a} 

This, as m the case of the corresponding toini in the longitude, 
is best considered in connexion with the elliptic inequality, and 
exactly the same results will follow 

Thus calling 3) , O, and a' the true longitudes of the moon, 
sun, and apse, the latter calculated on supposition of iinifoim 
motion, these two terms may he wntten, 

w = a[l f e cos( 3) “ a') + co3{ J) « a' -f 2(a' — o)]] 

= a[l 4 jE cos( T> -■ a' 4- S)] ; 
where -E cos S = e 4- ^me cos 2(a — O), 

-B sin 8 = sin 2{a' — 0)» 

These are identical with the equations of Art (70) 
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Vay lation 

86. To exjplain the effect of the term ra^a cos{(2 — 2m) 6 — 2^8}, 
w = a[l H- m“ cos{(2 — 2m) d — 2^8}] 

= a\l + rn? cos 2( J) — o)]. 

As far as this term is concerned^ the moon’s orbit would be 

an oval having its longest diameter m quadratures and least in 

syzygies Pnnaij^ia^ lib I , prop 66 cor 4 

The ratio of the axes of the oval oibit will be 

1 + . - , 

^^2 = §§ nearly, m being -0748. 

See Principia^ lib. Ill , prop 28 


Reduction 

87 The last important periodical term in the value of u is 

- — COS 2(5^0 -7). 

This, when increased by a constant^ is appioximately the dif- 
ference between the values of u in the orbit and m the ecliptic. 
For if be the reciprocal of the value of the radius vector 
m the orbit, 

Wj == w cos(latitude), 

“ 7(i+g^) = 

therefore w — mi^{gd — 7) 

= — \ali? cos 2 (^^ — 7) 

= const — cos2(;7^ — 7) 

88 The remaming terms m the value of u are of the thud 
ordei, and therefore very small* one of these corresponds to the 
annual equation m longitude Art (75), wheie it is of the second 
Older, having increased m the coiuse of mtegiation 
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Periodic time of the Moon 

89 We have seen, Ait (65), that the periodic tunc of tlie 
moon IS greater than if theio wcic no diRtmbuig toice , hut this 
refers to the mean pci iodic time cBtimatcd on an intcival of 
a great niunber of years, so that the ciioulai finietioim In the 
expression are then exticmcly small compaicd with the quantity 
which has nnifonnly inci cased 

When, however, wo considci only a few re\oluti(>ns, those 
terms may not all be neglected. The clhpttc 'ine<iuahft/ and 
the evection go through their values m aliout a month, the 
'oariation and reduction in about half-a-month , thoir elR'tts, 
theiefore, on the length of tific pciiod can seau‘cly be considered, 
as they will increase one poitiou and then deciease another of 
the same month 

But the annual equation takes one ycai to go thiough its 
cycle, and, duiing this time, the moon has disenbed ilinteen 
levolutions, hence, fluctuations may, and, as we sliall nowHlnwv, 
do take place in the lengths of the sidereal months duiuig the 
yeai. 

We have, considcimg only the annual ecjuation, Ait (75), 
j)t^Q sm {jnQ -f yS — f). 

Let The the length of the peiiod, thou when 6 is inm eased 
by 27r, t becomes t + 

thcrcfoic j)(^ + iZ^) = 27r -f 0 4* ?*iue hin(2>^w4- ntO ^), 
whence j? T = 277 4- Qune' am niir eos(>/i7r 4* 4- y8 — S’) ; 

theiefore T = mean period 4- sin }U7r eoH(o — f ), 

where O = mO 4- y8 4- niir = sun’s longitude at the bc'ginimig 
of the montli 4 mm- 

= aim’s longitude at the middle of tlu* month. 
Hence T will bo longest when o — 5* “ 9, 
and shortest when O — ^ = tt ; 

or jf’will be longest wlicn the sun at the middle of the month 
IS at perigee, and shoitest when at apogee ; hut, at pi ('sent, the 
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&*un IS at peiigee about the 1st of Januaiy, and apogee about 
the 1st of July; thciefoie, owing to annual equation^ the winter 
months will be longci than the summei months, the difference 
between a sideieal month in January and July, from this cause, 
being about 20 minutes 

90 All the inequalities or equations, which our expressions 
contain, have thus received a physical interpretation They 
wcie the only ones known befoie Newton had established his 
thcoiy^ but the necessity for such collections was fully iccog- 
nizod, and the values of the coefficients had already been pretty 
acciuatcly dctei mined; still, with the exception of the reduction, 
winch IS gcomctiically necessaiy, they were corrections em- 
piiically made, and it was scarcely to be expected that any 
but the larger inequalities, viz those of the first and second 
orders which wo have heie discussed, could be detected by 
obhoivation we find, however, that thieo otheis have, since 
Newton’s time, been indicated by obscivation before theory had 
<‘\plaxnod then cause These aic — the secular acceleration^ 
disco vexed by Halley, an inequality, found by Mayei, in the 
longitude of the moon, and of which the longitude of the ascend- 
ing node IS the aigumcnt, and finally an inequality discovered 
by Buig, which has only within the last five ycais obtained a 
solution For a further account of these, as also of some 
other inequalities which theoiy has made known, see Appendix, 
Alts (99), (100), (101), (102) 
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APPENDIX 

In this chapter will he found collected a few pioposltloiis 
mtliriately connected with the results oi the processes of the 
Lunar Theoiy as explained in the previous pages Reference 
has been made to some of them in the couise of the work, and 
the interest and importance of the othcis arc sufficient to justify 
their mtroduction here 

91 The moon is 7etained in hen orhit hy the foice of gravity^ 
that iis, iy the same force which acts on hodies at the siaface of 
the earth 

The proof of this is merely a numerical verification , the data 
required fiom observation are, 

the space fallen through from rest in l" by bodies at the oaith’s 
surface 16 1 feet, 

the radius of the earth ,, .. , . = 4000 miles, 

the periodic time of the moon . . == 27 J days, 

the distance of the moon from the earth’s cent! c= GO x 4000 miles. 

The force of the earth’s attraction cc Therefore, the 

, (dist )** ’ 

space fallen through in 1" at distance of moon by a body moving 

1 A 1 

from rest under the earth’s action = feet 

60 


= 00447 fret 
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But the moon m one second describes an angle 27 24 60**^ 
durmg which the approach to the earth 

= 60 X 4000 X 5280 (vers co) feet 
_ 60 X 4000 X 5280 27r^ . 

“ (27^)^ (24)^ (60r 

= *00448 feet 

TheieforCj the space through which the moon is deflected m 
one second from her straight path, is just the quantity thiough 
which she would fall towards the earth, supposmg her to be 
subject to the earth’s attraction, and we may, therefore, con- 
clude that she is retained in her orbit by the force of gravity 
When first Newton, m 1666, attempted to verify this result, 
he found a diffeience between the two values equal to one-sixth 
of the less , the reason of his failure was the mcorrect measures 
of the earth, which he made use of in his computation ; and it 
was not till about 16 years later that he was led to the true 
result, by usmg the more correct value of the earth’s radius 
obtained by Picart Pnncijpia^ lib iii , prop 4 


92 The moon^s orbit %s everywhere concave to the sun 

Let 8^ and M (fig 11) be the centres of the sun, earth, 
and moon We must bring the sun to rest by applying to the 
three bodies forces equal and opposite to those which act on 
the sun ; but these are so small that we may neglect them and 
consider the moon as moving round the sun fixed, and disturbed 
by the earth alone 

The forces on M are, therefore, m M8^ 

E 


and 




mME 


This last must be resolved into two, one in -MS, the other pei- 
pendicular to it 

Therefore, the whole central force on the moon in 

77l' E 
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and the proposition will be proved if we show that this force is 
always positive 


27r SG~^ 8]\r 


Now, period lound snn = — ^ — = 


m ■ 


neaily, 


and 

therefoie 

therefore 


therefore 


earth = 


27r Al/" 


m' E + M , 

nearly, 

E 


m 


> 


SM^ " EM'^ ’ 

m' 8M E 

8M^ EM EM‘ 


m 

'm 


E 




E 

ME ’ 


but the least value of the ccntial foicc coiicsponds to cosil/= — I, 

and is then It is, thciofoic, always positive, or 

the path always concave to the sun 

At new moon the foico with winch the moon tennis to the 
sun IS, theiefoie, greater than that with whu‘h she tennis to tlic 
earth the earth being itself m motion in the same ilnoction, 
and, at that instant, with gieatci velocity, will easily explain 
how, notwithstanding this, the moon still revolves about it. 


Cential and Tangential DiBtuihing Fotees 

93 We have hithei’to considoicd the effe(‘ts of the eentral 
and tangential distuibing foiccs in combination; but it will 
be interesting to determine to which of them tlic siweial in- 
equalities pnncipally owe their existence 

(1) To dete'^mine the effect of the central disfurhnu/ force. 

Make T = 0 ; 

.1 , d^if P 

theiotoie ^ 0 

dd Iru" ^ 
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and tins being substituted in the diffeiential equation (B') of 
Art (30), gives 

d‘u 

.. . + Acos(0-a)+ .. , 

tlie solution of which is 

= a{l + <3cos((9-.a)} + ... -f 

Our fiist approximate value w = a {1 + ecos(0~a)} is thus 
corrected by a teim which, on account of the factor 0, admits 
of indefinite increase, and thus becomes ultimately a more 
impoitant teim than that with which we started as bemg very 
nearly the true value, and which is confirmed as such by obser- 
vation (22) : for the moon’s distance, as determmed by her 
paiallax, is never much less than 60 times the earth’s radius ; 
whereas this new value of w, when 6 is very great, would make 
the distance indefinitely small, and, on the same principle, 
we see that any solution, which comprises a term of the form 
A6 sin(0 — a), cannot be an approximate solution except foi 
a small range of values of 0. 

Such terms ‘if they really had an existence in our system, must end 
‘ in its destruction, or at least m the total subversion of its present state , 

‘ but when they do occur, they have their origin, not m the nature of the 
‘ differential equations, hut m the imperfection of our analysis, and m the 
‘ inadequate lepresentation of the perturbations, and are to he got rid of, 

‘ or rather included in more general expressions of a periodical nature, by 
‘ a more refined investigation than that which led us to them The nature 
‘ of this difficulty will be easily understood from the foUowmg reasoning 
‘ Suppose that a tem, such as a sin {AO f JB), should exist m the value 
‘ of in which A being extremely minute, the period of the inequality 
‘ denoted by it would be of great length, then, whatever might be the value 
‘ of the coefficient a, the inequality would still be always confined within 
* certain limits, and after many ages would return to its former state 

‘ Suppose now that our peculiar mode of arriving at the value of w, led 
‘ us to this term, not in its real analytical form a sin {AO + B)^ hut by the 


IS the equation for determiiimg the constant e, and m the second approxima- 
tion, « would be found from 

0=s« + 2fi sm(« — a) + fe® sra2 («— a)-j- , 

giving different values of e at each successive approximation 
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* waj of Its development in powers ot 0, a \ ^0 { 7 ^ 1 A c , and that, not 
‘ at once, hut piecemeal, as it were, a fust approximation giving us onl} 
‘ the term a, a second adding the term /30, and so on If w e stopped here, 
‘ it IS ob\ious that we should mistake the nature of this mequahty, and that 
‘ a really periodical function^ from the effect of an imperfect approximation* 

* would appear under the form of one not periodical These terms 

* m the value of u, when they occur, a^’e not superfluous , they are essential 
‘to its expression, but they lead us to eironeous conclusions as to the 
‘ stability of our system and the general laws of its pertuibations, unless 
‘ we keep in view that they are only pat ts of senes , the principal parts, 
‘ it is true, when we confine ourselves to intervals of moderate length, 
‘ but which cease to be so after the lapse of very long times, the rest of 
‘ the series acquiring ultimatelj the preponderance, and compensating the 
‘want of periodicity of its first terms’ — Sir John Hleschel, Ency- 
clopcedia Mettopoktana — Physical Astronomy, p 679 


36 To extricate ourselves from this difficulty,' and to alter 
the solution so that none but periodical teims may be intro- 
duced, let us again observe that the equation ^ -f w = = a 


of Art (30), which gave the solution and thus led to the 
difficulty, IS only an approximate form of the fiist order of the 
exact equation (^S') of the same article. Any value of w, therefoie, 

d^u 

which satisfies the approximate equation ^ -h = a to the 


first order, and which evades the difficulty mentioned above, 
may be taken as a solution to the same order of the exact 
equation (^'). 


Such a value will be 


= a {1 + e cos (c0 — a)} ?7/, 

piovided 1 — o'* be of the first order at least, for then 
d^u 

— + w = a + ae (1 — o'*) cos [cd — a), 

= a to the first order 


37. The observations hitherto made to check our approx- 
imations weie extremely rough (22), and carried on only for 
a short interval; but when they aie made with a little more 
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accuiacy, and extended over seveial revolutions of the moon, 
it IS found that her apse and the plane of hei orbit aie m 
constant motion 

The above form of the value of u is suggested by our pie- 
vious knowledge of this motion of the apse, which, as we shall 
see Art (66), is connected with the value of c here introduced; 
and there is no doubt that Clairaut, to whom this artifice is 
due, was led to it by that consideration, and by his acquaint- 
ance with the results of Newton’s ninth section, which, when 
translated into analytical language, lead at once to this foim 
of the value of w* 

We might, theiefore, taking for granted the results of ob- 
servation, have commenced our approximations at this step, 
and have at once written down 


u = a {1 4“ e cos {c6 - a)}, 


but we should, m so doing, have merely postponed the difficulty 
to the next step, since there again, as we shall find, the dit- 
ferential equation is of the foim 


d^u 

W 


4- w = a function of 6, 


the collect integial of which would be, 

w = ^ cos (0 — jB) H- , 

and this would at the next opciation bring in a teim with 9 foi 
a coefficient, which we now know must not be We shall, therc- 


* Newton has there shewn, that if the angular velocity of the oihit he 
to that of the body as G - F to G, the additional centiipetal force is 
- F^ 

— — 4V, the oiigmal foice being jllu* Therefore 


d^u fji 

A* 


GP-~F^ 


w, 


dif G‘ “ W ~ Q‘ h,‘ CP 


a e cos 




where is the same as oui c 


D 




